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OPERADORES ORTOGONAIS 


Ao final desta aula, você deverá ser capaz de: 


compreender o conceito e as propriedades apresentadas 
sobre operadores ortogonais; 

aplicar os conceitos apresentados em exemplos impor- 
tantes. 


Álgebra Linear II | Operadores Ortogonais 


Pré-requisitos 


Aulas 10 a 14,17 e 


18. 
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OPERADORES ORTOGONAIS 


Você deve se lembrar de que um operador linear T :R” — R” é dito 
ortogonal se existe uma base ortonormal & de IR” tal que a matriz de T 
na base & é uma matriz ortogonal, isto é, se a matriz [T|w é ortogonal. 


Veremos que os operadores ortogonais estão bem definidos no sen- 
tido de que o fato de ser um operador ortogonal não depende da base 
ortonormal escolhida, ou seja, se a matriz [T |, numa certa base ortonor- 
mal q de R”, for ortogonal, então a matriz [7] B também será ortogonal 
para qualquer outra base ortonormal B de R”. 


Na verdade, temos o seguinte resultado: 


Teorema 19.1. 


Sejam T :R” — R” um operador linear ortogonal e a e B duas bases 
ortonormais de R”. Se a matriz [T]w é ortogonal, então a matriz [T]p 
também será ortogonal. 


Demonstração 


O teorema sobre mudança de base para operadores lineares, visto no 
curso de Algebra Linear I, nos garante que 


[Tp EA 


onde P é a matriz mudança de base, da base B para a base x. Como 
o e B são duas bases ortonormais de R”, temos que P é uma matriz 
ortogonal e, pelo Teorema 10.1. da Aula 10, segue-se que 


pPl=Pp', 
onde P' é a transposta da matriz P. Assim, 


[7]; = PT)oP. 


Como [T|w é uma matriz ortogonal por hipótese e como o produto 
de matrizes ortogonais é também uma matriz ortogonal, concluímos que 
[7] também será uma matriz ortogonal. 


O resultado anterior simplifica um problema crucial: para verificar- 
mos se um dado operador linear T: R” — Rºé ortogonal, basta consi- 
derar qualquer base ortonormal a de R” e verificar se a matriz [T|o é 
uma matriz ortogonal. 


Exemplo 19.1. 


Verifique que o operador linear T:Rº* > Rº 


T(x,y,2) = (xcos6 — ysen O, xsen6 +ycos6, 2), 


com 0 € [0,27), é um operador ortogonal. 


AULA ES MÓDULO 2 


Solução: 
De fato, escolhendo a base canônica fe,,e>,e3! de Rº, dada por 
e/=(1,0,0), e-=(0,1,0) ces = (0,0,1), 
obtemos 
T(ei) = (cos 0,sen 0,0) 


T(e>) = (— sen8 cos0,0) 
T(es) = (0,0,1). 


Portanto, a matriz que representa T nesta base é dada por 


cos0 -—senO O 
A= sen O cos0 O 
0 0 1 


Sabemos que A é uma matriz ortogonal de Rê. Mais ainda, A é uma rotação 
de O radianos em torno do eixo-z (Exemplo 17.1 da Aula 17). Assim, o opera- 
dor linear T é um operador ortogonal. 


O próximo teorema segue imediatamente do Teorema 10.2 da Aula 
10. 


Teorema 19.2. 


Seja T:R” — R” um operador ortogonal. Então as seguintes pro- 
priedades são válidas: 


1. T transforma bases ortonormais em bases ortonormais, ou seja, se 
(Vi,V2,...,Vny é uma base ortonormal de IR”, então 
(Tvi,Tvo,...,Tvn+ também é uma base ortonormal de R”. 


2. T preserva o produto interno, ou seja, para todo u,v € R” vale 
que 


(Tu,Tv) = (u,v). 
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3. T preserva a norma, ou seja, para todo v € R” vale que 


EP vil = Ihvll. 


Exemplo 19.2. 


Seja T : IR? — R? um operador ortogonal, então sua matriz na base 
canônica é da forma 


cos0 -—sen6 a cos 0 sen 0 
sen 6 cos 6 E senO -—cos0 /' 


onde 6 E [0, 27). 


Solução: 


De fato, sendo T : R? — R? um operador ortogonal, sua matriz na base 
canônica de R? será uma matriz ortogonal de ordem 2. Mas, pelos Exemplos 
10.1 e 10.2 da Aula 10, sabemos que toda matriz ortogonal de ordem 2 é da 


forma 
cos0 —senô cos 9 sen 6 
ou 
sen 6 cos 6 sen6 —cos0 
Sabemos também que a primeira matriz representa uma rotação de O radi- 
anos, no sentido anti-horário, em torno da origem, e a segunda matriz repre- 


senta uma reflexão em torno da reta pela origem que forma um ângulo de 6/2 
radianos com o semieixo x positivo. 


Exemplo 19.3. 


a. Determine a transformação linear T :R? — Rº que leva o seg- 
mento de reta de extremidades (—6,2) e (—1,2) ao segmento 
de reta de extremidades (—2,6) e (1,2), respectivamente (veja 
a Figura 19.1). 


b. Mostre que a transformação acima é uma rotação. Determine, 
também, o ângulo dessa rotação. 


AULA ES MÓDULO 2 


Figura 19.1: O operador 7. 
Solução: 


a. Queremos encontrar escalares a,b,c,d E R tais que a matriz que repre- 
senta T na base canônica seja dada por 


m=(18) 


Da condição sobre as extremidades, temos 


É fácil ver que a solução desse sistema é dada por: 


a=3/5; b=4/5; c=-4/5 e d=3/5. 


Assim, 
o 3/5 4/5 
E] = -4]s 3/5 )- 


b. Como as colunas da matriz [7], representadas pelos vetores 
vi=(3/5,-4/5) e v> = (4/5, 3/5), formam uma base ortonormal de 
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R?, concluímos que a matriz [7] é ortogonal e, consequentemente, o 
operador linear T é um operador ortogonal. Além disso, det/T] = 1 e, 
assim, o operador T é uma rotação de R? cujo ângulo O é dado por 


0 = —arccos(3/5). 


Exercício 19.1. 


I. Seja T:Rº — Rº uma reflexão num plano x de Rº tal que 
T(1,0,-1) = (—1,0,1). Determine a matriz que representa o 
operador T com respeito à base canônica. 


2. Determine os autovalores e os autovetores associados da transfor- 
mação linear T do exercício anterior. 
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Aula 20 : 


PROJEÇÕES ORTOGONAIS — 1º PARTE 


Ao final desta aula, você deverá ser capaz de: 


compreender o conceito de projeção ortogonal em di- 
mensão 2; 

aplicar os conceitos apresentados em exemplos impor- 
tantes. 


Álgebra Linear II | Projeções ortogonais — 1º Parte 


Pré-requisitos 
Aulas 10 a 14, 17, 
18e 19. 
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PROJEÇÕES ORTOGONAIS — 1º PARTE 


Nesta e na próxima aula, vamos apresentar um tipo de transforma- 
ção usada em áreas como a Computação Gráfica e o Desenho Geométrico. 
Trata-se das projeções ortogonais. Nesta primeira aula, trabalharemos 
com as projeções ortogonais em R2. 


Exemplo 20.1. 


Determine a matriz que representa a projeção ortogonal sobre o 
eixo-x, isto é, sobre a reta de equação cartesiana y = 0. 


Solução: 


Geometricamente, essa transformação é representada pela Figura 20.1. 


T(w=(x,0) 


Figura 20.1: A projeção ortogonal no eixo-x. 


Assim, temos a transformação linear 


T:R?2SR? 
T(x,y) = (2,0). 


Denotando por fe,,e>) a base canônica de R?, temos que 
T(e)=T(L,0)=(1,0)=1-e,+0-e, 


T(e,) = T(0,1) — (0,0) =0-e;+0-e. 


Portanto, a matriz que representa a transformação T na base canônica é 


1=(00) 


dada por 


Vemos imediatamente algumas propriedades dessa projeção ortogonal. 


1. A matriz A e, portanto, o operador T, não é inversível, pois det(A) = 0. 


2. Como T(e,;) = 0- e», então À» = O é um autovalor de T com autovetor 
associado e; = (0,1). Não é difícil ver que o autoespaço associado a 
À» = 0 é exatamente o eixo-y, isto é, a reta de equação cartesiana x = 0. 


AULA ES MÓDULO 2 


3. Como T(e;) = 1-e;, então À = 1 é um autovalor de T com autovetor 
associado e; = (1,0). Não é difícil ver que o autoespaço associado a 
AM = 1 é exatamente o eixo-x, isto é, a reta de equação cartesiana y = 0. 


4. O operador T é diagonalizável e seu polinômio característico é 
p(x) =x(x— 1). 


Exemplo 20.2. 


Determine a matriz que representa a projeção ortogonal sobre o 
eixo-y, isto é, sobre a reta de equação cartesiana x = 0. 


Solução: 
A projeção ortogonal no o eixo-y é dada pela transformação linear 
T:RSR? 
T(x,)) = (0,9). 


Geometricamente, esta transformação é representada pela Figura 20.2. 


Tisy)dl.. 


Figura 20.2: A projeção ortogonal no eixo-y. 
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Como no Exemplo 20.1, temos que 


nn 
PE 
asa 
— 
I 
a 
A 
a 
[= 
— 
I 


(0,0) =0-e/+0-e> 
(0,1)=0-ej+1-e,. 


a 
PEN 
o 
2) 
— 
I 
a 
EN 
e) 
ta 
— 
I 


Portanto, a matriz que representa a transformação T na base canônica é 
dada por 


s=(8 8) 


Como antes, vemos que: 


« A matriz A e, portanto, o operador T, não é inversível, pois det(A) = 0. 


. Como T(e;) = 0-e;, então À = O é um autovalor de T com autovetor 


associado e; = (1, 0). Não é difícil ver que o autoespaço associado a 
AM = O é exatamente o eixo-x, isto é, a reta de equação cartesiana y = 0. 


. Como T(e,) = 1 -e>, então À, = 1 é um autovalor de T com autovetor 


associado e; = (0, 1). Não é difícil ver que o autoespaço associado a 
À» = 1 é exatamente o eixo-y, isto é, a reta de equação cartesiana x = 0. 


2 


4. O operador T é diagonalizável com polinômio característico 


p(x) =x(x— 1). 


Os Exemplos 20.1 e 20.2 são muito simples, porém são muito im- 
portantes a sua compreensão e o seu significado geométrico. Especial- 
mente, certifique-se de que tenha entendido os autoespaços associados 
a cada autovalor. Usaremos essas ideias para apresentar a projeção or- 
togonal sobre uma reta L qualquer do R? passando pela origem. Se você 
compreendeu bem a geometria dos exemplos anteriores, então não terá 


dificuldade em acompanhar o caso geral a seguir. 


Exemplo 20.3. 


Descreva a projeção ortogonal sobre uma reta L de 
pela origem. 


Solução: 


R? que passa 


Suponhamos que a reta L seja paralela a um vetor unitário u, € R?, como 


ilustra a Figura 20.3. 


Figura 20.3: A reta L paralela ao vetor unitário u,. 


O efeito geométrico da projeção ortogonal sobre a reta L é observado na 
Figura 20.4. 


Figura 20.4: A projeção ortogonal na reta L. 


A projeção ortogonal de um vetor v na direção do vetor u; é dada por 


TR? > R 
(v,u1) 


Vs Dai u,, 


de onde vemos que T é uma transformação linear. Para obter a fórmula acima, 
observamos que desejamos um vetor Tv da forma Tv = ku, de modo que 
v— ku; seja ortogonal au, , como indica a Figura 20.5. 


AULA ES MÓDULO 2 
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Figura 20.5: A projeção ortogonal de v na direção de u; 


Assim, da ortogonalidade entre vy — ku, e u:, temos 


O = (v-ku,,u;) 
= (v,u1) — (ku,,ur) 
= (vu))—k(u,,u)), 


o que nos dá 


k(u,,u,)) = (vu) 
k = (v,u1) : 
(u,,U1) 
e, portanto, 


(1,41) 


Observe que na fórmula acima o vetor u; não precisa ser unitário, mas, 
caso seja, como (u,, u1) = 1, então a fórmula acima se simplifica para 


Tv=(v,u;) ui. 


Nosso problema agora é encontrar a matriz que represente a transformação 
T. Veremos que, escolhendo uma base ortonormal adequada de RZ, a matriz de 
T nessa base é muito similar à matriz do Exemplo 20.1, visto anteriormente. 
Lembre que o problema da escolha de uma base ortonormal adequada já foi 
tratado quando estudamos as reflexões de R? com respeito a uma reta qualquer 
passando pela origem. Veja a Aula 12. 


Seja 8 = fu,,u>) uma base ortonormal de R? onde u, é um vetor unitário 
paralelo à reta L e u> é um vetor unitário normal à reta L. Veja a Figura 20.6. 


Nesse caso, como (u,,u,) = 1 e pela observação acima, temos que 
Tv = (v,u/) u;. Assim, vemos que 


Tu =(u,u))u=u=1-u,+0-u> 
Tu = (u,u,;) u, =0-u, =0-u, +0-u,». 


Figura 20.6: A base ortonormal 8 = (u,, u2>. 


Portanto, a matriz que representa a transformação T na base B é dada por 


r=(0 0): 


que é exatamente da mesma forma que a matriz do Exemplo 20.1. Se qui- 
sermos obter a matriz que representa T na base canônica, é só fazermos uma 
mudança de base. Se a = fe1,e>) é a base canônica de Rº, então 


[Ta E P[T]gP”!, 


onde P é a matriz mudança de base, da base B para a base q. Como 
P=[u, us), isto é, suas colunas são vetores ortonormais, então P é uma matriz 
ortogonal e, portanto, P-! = P'. Como nos Exemplos 20.1 e 20.2, temos as 
seguintes propriedades. 


1. As matrizes [7]q e [T]p e, portanto, o operador T, não são inversíveis, 
pois det [7]p = 0. 


2. Como T(u,) = 0-u», então À, = O é um autovalor de T com autovetor 
associado u2. Não é difícil ver que o autoespaço associado a À, = O é 
exatamente a reta pela origem ortogonal à reta L. 


3. Como T(u,) = 1-u;, então À, = 1 é um autovalor de T com autovetor 
associado u,. Não é difícil ver que o autoespaço associado a À = 1 é 
exatamente a reta L. 


4. O operador T é diagonalizável e seu polinômio característico é 
plo) =x(x— 1). 


Cabe aqui, mais uma vez, ressaltar a analogia entre este terceiro exemplo 


AULA ES MÓDULO 2 


CEDERJ 19 


Álgebra Linear II | Projeções ortogonais — 1º Parte 


e os dois primeiros. Isto se deve à escolha adequada de uma base ortonormal 
de R2. 


Exercício 20.1. 


1. Determine a matriz da projeção ortogonal sobre a reta y = 3x 
com respeito à base canônica. 


2. Determine os autovalores e os autoespaços associados da trans- 
formação linear do exercício anterior. 
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Aula ) | 


PROJEÇÕES ORTOGONAIS — 2º PARTE 


Ao final desta aula, você deverá ser capaz de: 


compreender o conceito de projeção ortogonal em di- 
mensão 3; 

aplicar os conceitos apresentados em exemplos impor- 
tantes. 
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PROJEÇÕES ORTOGONAIS — 2º PARTE 


Pré-requisitos 

Aulas 10a 14,17a Nesta aula daremos continuidade ao estudo das projeções ortogo- 

20. nais, estudando as projeções ortogonais em R?. Apresentamos, ini- 
cialmente, os casos mais simples das projeções ortogonais nos planos 
coordenados. Em seguida, trataremos do caso geral de uma projeção 
ortogonal sobre um plano passando pela origem. 


Exemplo 21.1. 


Determine a matriz que representa a projeção ortogonal sobre o 
plano-xy, isto é, sobre o plano de equação cartesiana z = 0. 


Solução: 


Geometricamente, essa transformação é representada pela Figura 21.1. 


V= (6360) 


» 


ve (x,,0) 


X 


Figura 21.1: A projeção ortogonal no plano-xy. 


Assim, temos a transformação linear 


T:R SR? 
T(x,7,2) = (x,7,0). 


Denotando por fe;,e>,e3! a base canônica de Rº, temos que 


T(e1) =T(1,0,0) = (1,0,0) = 1 .e/+0-e2+0-e3 
T(e)=T(0,1,0)=(0,1,0)=0-e;+1-e2+0-e; 
T(es) = T(0,0,1) = (0,0,0) =0-e,+0-e+0-es. 


Portanto, a matriz que representa a transformação T na base canônica é 
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dada por 
100 
A=|0 10 
000 


Como nos exemplos da Aula 20, vemos imediatamente algumas propriedades 
dessa projeção ortogonal. 


AULA E MÓDULO 2 


1. A matriz A e, portanto, o operador T, não são inversíveis, pois 
det(A) = 0. 


2. Como T(e3) = 0-es, então Às = O é um autovalor de T com autovetor 
associado e3. Não é difícil ver que o autoespaço associado a Às = O é 
exatamente o eixo-z, que é o espaço gerado por es. 


3. ComoT(e;))=1-ejeT(e,)=1-e,, então À = À» = 1 é um autovalor 
de T de multiplicidade 2 com autovetores associados e; e e,. Não é 
difícil ver que o autoespaço associado a À, = À» = 1 é exatamente o 
plano-xy, que é o espaço gerado pelos vetores canônicos e; e ep. 


4. O operador T é diagonalizável com polinômio característico 
pl) =x(x—1)2, 


Mais uma vez, chamamos a atenção do aluno para que compreenda 
bem a geometria desse exemplo, pois ela será recorrente nos exemplos 
seguintes. Vejamos outro exemplo de projeção ortogonal em um plano 
coordenado. 


Exemplo 21.2. 


Determine a matriz que representa a projeção ortogonal sobre o 
plano-yz, isto é, sobre o plano de equação cartesiana x = 0. 


Solução: 


Geometricamente, essa transformação é representada pela Figura 21.2. 
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PRERCRRRI gare ee  V=(0,7,2) 


X 


Figura 21.2: A projeção ortogonal no plano-yz. 


Assim, temos a transformação linear 


T:RS Rê 
T(x,7,2) = (0,9,2). 


Se você entendeu bem a geometria do Exemplo 21.1, então verá que, neste 
caso, temos 


T(e1) = T(1,0,0) = (0,0,0) =0-e; +0-e, +0-es 
T(e») = T(0,1,0) == (0,1,0) =0.e;+1 -e+0-e3 
T(es) = T(0,0,1) == (0,0,1) =0-e; +0-e> +1 -es. 


Portanto, a matriz que representa a transformação T na base canônica é 
dada por 


O 00 
A=| 010 
001 
Seguem também as propriedades: 


1. A matriz 4 e, portanto, o operador T, não são inversíveis, pois 
det(A) = 0. 


2. Como T(e;) = 0-e;, então À, = O é um autovalor de T com autovetor 
associado ej. Não é difícil ver que o autoespaço associado a À = O é 
exatamente o eixo-x, que é o espaço gerado por e. 


3. Como T(e;)=1.ezeT(es)=1-es, então À = Às = 1 é um autovalor 
de T de multiplicidade 2 com autovetores associados e, e es. Não é 
difícil ver que o autoespaço associado a À, = À3 = 1 é exatamente o 
plano-yz, que é o espaço gerado pelos vetores canônicos e, e es. 
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4. O operador T é diagonalizável com polinômio característico 
p(x) =x(x— 1). 


O outro caso trivial, a projeção ortogonal sobre o plano-xz, é total- 
mente análogo aos exemplos anteriores e deixamos como exercício para 
você. Assim, estando bem compreendidos os dois exemplos anteriores, 
podemos tratar da projeção ortogonal sobre um plano qualquer de Rº? 
passando pela origem. 


Exemplo 21.3. 


Descreva a projeção ortogonal sobre um plano x de Rº que passa 
pela origem. 


Solução: 


Seja T : Rê — Rº a projeção ortogonal sobre o plano x. Geometricamente, 
essa transformação é representada pela Figura 21.3. 


Plano 7 


Figura 21.3: A projeção ortogonal no plano-7. 


Vamos agora obter uma base ortonormal 8 de IR? de modo que a matriz que 
representa a transformação T nessa base seja da mesma forma que a matriz do 
Exemplo 21.1. Como conhecemos a equação cartesiana de plano 7, sabemos 
como obter um vetor normal a esse plano. Lembre: se x tem equação ax +by + 
cz+d =0, então o vetor u = (a,b,c) é um vetor normal ao plano 7. Seja, então, 
us um vetor unitário normal ao plano 7. Usando a equação cartesiana de 7, 
como foi feito nas Aulas 17 e 18, facilmente determinamos vetores unitários uU; 
e u» de modo que 8 = fu,,u>,u3) seja uma base ortonormal de R?. Observe 
que os vetores unitários u; e u; são ortogonais e pertencem ao plano 7. 


Veja a Figura 21.4. 
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Us 


Uh Plano Z 


Figura 21.4: A base ortonormal 8 = (u,,u>,us). 


A projeção ortogonal de um vetor v sobre o plano 7 é dada por 


TR Ss R 
(v,u1) (v,u2) 


vo Tv= 1 
(u,,U1) (u>,U2) 


Us, 


de onde vemos que T é uma transformação linear. Para obter a fórmula acima, 
observamos que desejamos um vetor Tv da forma Tv = ku; + kou> de modo 
que v— ku, — ku; seja ortogonal a u; e us, como indica a Figura 21.5. 


lano x 
V -ku,- ku, P 
Ug2 
v 
kju+Ius 
us 


Figura 21.5: A projeção ortogonal de v no plano 7. 


Assim, da ortogonalidade entre vy — ku; — ku, e us, temos 


= 


v—ku; — ku,,u)) 
v,U1) id (ku, u,) = (kouo,U|) 
v, u,) — ki (u,,U1) — ko (u>,U1) 
Vv, u,) — ki (u,,U1) , 


( 
( 
( 
( 


já que (u5,u,) = 0, o que nos dá 
ki (1,41) = (v,u1) 


E (v,U1) 
(ui). 
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Analogamente, da ortogonalidade entre v — ku, — ku, e us obtemos que 


Vo) 
(uo,u2)” 


e, portanto, 


v,u v,u 
Tv=ku,+hou = po Gis; A) u 
(u,,U1) (u2,U2) 
Usando o fato de u, e u; serem vetores unitários, isto é, (u,,u,) = (U5,U2) = 
1, obtemos 
Tv=(v,u,)u;+(v,u2) us. 


Portanto, vemos que 


Tu, = (u,,u,) u+(u,,u2) w =u,=1-u,+0-u2+0-us 
Tu, = (u,,u,) u, + (u2,u2) uw=u=0-u,4+1-u5+0-us 
Tu; = (us,u,) uj + (us,u2) u2 — 0 -u; +0-u2+0-us. 


Portanto, a matriz que representa a transformação T na base B é dada por 


que é exatamente da mesma forma que a matriz do Exemplo 21.1. Se qui- 
sermos obter a matriz que representa T na base canônica, é só fazermos uma 
mudança de base. Se « = [e,,e>,€3) é a base canônica de RR, então 


Fla=PIrpr 


onde P é a matriz mudança de base, da base B para a base q. Como 
P=[u, u2 us], isto é, suas colunas são vetores ortonormais, então P é uma 
matriz ortogonal e, portanto, P”! = P”. Como nos exemplos 21.1 e 21.2, temos 
as seguintes propriedades: 


1. As matrizes [T]y e [7], e, portanto, o operador T, não são inversíveis, 
pois det [7]p = 0. 


2. Como T (us) = O-us, então Às = O é um autovalor de T com autovetor 
associado u3. Não é difícil ver que o autoespaço associado a Às = O é 
exatamente a reta pela origem ortogonal a 7. 


3. Como T(u;))=1-u/ e T(u,) = 1-u,, então À = À» = 1 é um auto- 
valor de T com autovetores associados u, e us. Não é difícil ver que o 
autoespaço associado a À, = À, = 1 é exatamente o plano 7. 
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Z 


4. O operador T é diagonalizável com polinômio característico 
po ta DA 


Cabe aqui, mais uma vez, ressaltar a analogia entre este terceiro 
exemplo e os dois primeiros. Isso se deve à escolha adequada de uma 
base ortonormal de Rº. 


Exercício 21.1. 


1. Determine a matriz da projeção ortogonal sobre o plano-xz com 
respeito à base canônica. 


2. Determine a matriz da projeção ortogonal sobre o plano x — z = 0 
com respeito à base canônica. 


3. Determine a matriz da projeção ortogonal sobre o plano gerado 
pelos vetores v; = (1,1,0) e v2 = (—1,1,1), com respeito à base 
canônica. 
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Aula ) ) 


MATRIZES SIMÉTRICAS 


Ao final desta aula, você deverá ser capaz de: 


compreender o conceito de matriz simétrica; 


aplicar os conceitos apresentados em exemplos impor- 
tantes. 


Álgebra Linear II | Matrizes Simétricas 


Pré-requisitos 
Aulas 6, 7, 8,9, 10, 
20€e21 
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MATRIZES SIMÉTRICAS 


Em muitas aplicações da Álgebra Linear, as matrizes simétricas 
aparecem com maior frequência que qualquer outra classe de matrizes 
importantes. A teoria correspondente a essas matrizes é muito rica e 
elegante, e depende, de maneira especial, das teorias de diagonalização 
e ortogonalidade, vistas em aulas anteriores. Veremos, nesta aula, que 
a diagonalização de uma matriz simétrica é um fundamento essencial 
e necessário à discussão das formas quadráticas que estudaremos no 
próximo módulo. 


Lembramos que todas as matrizes e vetores considerados têm so- 
mente elementos e componentes reais. Antes de começarmos a estudar 
a teoria de diagonalização de matrizes simétricas, convém lembrarmos 
de algumas definições que serão essenciais a este conteúdo. 


Definição 22.1. 


Uma matriz A E Ma(R) é simétrica se A = A, onde A” representa 
a matriz transposta de 4. Equivalentemente, a matriz A = (ai;) é 
simétrica se a;; = a j; para todo i, j. 


Observe, primeiramente, que o conceito de matriz simétrica se aplica 
apenas a matrizes quadradas. Observe também que os elementos da 
diagonal principal de uma matriz simétrica A podem assumir valores 
arbitrários; no entanto, elementos simétricos com respeito à diagonal 
principal têm o mesmo valor. 


Exemplo 22.1. 


As duas matrizes a seguir são simétricas: 


4-1. 0 
a=(13) e B=|.=1 2 3 
Rj! da 2 


No entanto, as matrizes abaixo não são simétricas: 


= qdo s= 
0] e D= 42 
e 


A matriz C não é simétrica porque ela não é matriz quadrada, e a 
matriz D não é simétrica porque dy, = 14 —1 = dis. 
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Vamos rever algumas propriedades das matrizes simétricas. 
Teorema 22.1. 


Sejam 4,B E M,(IR) matrizes simétricas. Então A4+B e c4, onde 
c E IR, também são matrizes simétricas. 


Vale observar que o produto de duas matrizes simétricas não é ne- 
cessariamente uma matriz simétrica. Por exemplo, dadas as matrizes 


simétricas 
1: :2 45 
ie Re 
temos que a matriz produto 
1 3 45 I4 17 
aB= (3 :) e elo Sl 
não é uma matriz simétrica, pois (AB)21 = 23 £ 17 = (AB) 2. 


Vamos rever o processo de diagonalização de matrizes, descrito 
nas Aulas 6 e 7, agora aplicado a um caso particular de uma matriz 
simétrica. 


Exemplo 22.2. 


6 —2 —1 
Diagonalize, caso seja possível, amatrizA= | —2 6 — 
-1 -1 5 
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Solução: 


O polinômio característico da matriz A é dado por: 


p(x) = det(x—A) 
x— 6 2 1 
= 2 x—6 1 
1 1 x—s5 
x— 6 1 2 1 2. 1 
= fgi6): 1 x—s5 [2 1 x—S5 ESA x—6 1 | 


22 — 17x? +90x— 144. 


As possíveis raízes racionais de p(x) são, obrigatoriamente, divisores de 
144. Por inspeção, vemos que 3 é uma raiz e, depois, completando a fatoração 
de p(x), descobrimos que 6 e 8 também são raízes. Assim, 


pl) = (x—3)(x— 6) (x — 8). 


Assim, os autovalores da matriz A são A — 3, App — 6 e À — 8. Como 
a matriz A possui 3 autovalores distintos, já podemos concluir que ela é uma 
matriz diagonalizável. 


Para o autovalor A, = 3, temos que os seus autovetores associados, 
v= (x,,2), satisfazem o sistema linear 


(35 — A)v = 0. 


Um cálculo rotineiro, como foi visto na Aula 7, mostra que o autoespaço 
E(3) é um subespaço de dimensão 1 e é gerado pelo vetor v; = (1, 1, 1). Ana- 
logamente, o autoespaço E (6), associado ao autovalor À» = 6, é o subespaço 
de dimensão 1 gerado pelo vetor vz = (—1,-—1,2), e o autoespaço E (8), as- 
sociado ao autovalor Às = 8, é o subespaço de dimensão 1 gerado pelo vetor 
v3 = (—1,1,0). Esses três vetores, v1, V> e v3, formam uma base de Rº e pode- 
riam ser usados para construir uma matriz P que diagonaliza a matriz A. É 
fácil ver que fvi,v2,V3) é um conjunto ortogonal de Rº e que obteremos uma 
matriz ortogonal P se usarmos uma base ortonormal fu,,u2,us), obtida de 
(vi,Y2,V3 |, normalizando cada um dos vetores v1, Y2 e v3. Como um múltiplo 
não-nulo de um autovetor também é um autovetor, a nova base (u,,u2,us) 
também seria uma base de autovetores de IRÊ. Os vetores assim obtidos são: 


u=(1/v3,1/V3,1/v3); 
w=(—1//6,-1//6,2/V6) e 
us=(—1/v2,1//2,0). 


Assim, as matrizes P e D são dadas por: 


1/43 —1//6 —1/v2 
P=| 1/03 -1//6 1/42 e D= 
1/3 2//6 0 


[em a SO) 
Sano 
o So o 


Sabemos, das Aulas 6 e 7, que A = PDP”. Agora, como as colunas de 
P formam vetores ortonormais, então, pelo Teorema 9.2 da Aula 9, P é uma 
matriz ortogonal, isto é, P! — P'. Assim, temos também que A = PDP". 


Vimos, no Exemplo 22.2, que os autovetores da matriz simétrica 
A, associados a autovalores distintos, são ortogonais. Isso é uma pro- 
priedade geral, como mostra o próximo teorema. 


Teorema 22.2. 


Seja A E M,(R) uma matriz simétrica; então qualquer conjunto de 
autovetores associados a autovalores distintos são ortogonais. 


Demonstração 


Sejam Vi, V2,...,v; autovetores da matriz A associados aos autova- 
lores distintos A1,Ã»,...,Ax. Assim, dados À; 4 À;, e observando que 
Av;=Ajv;e Av; = AN js queremos mostrar que (Vi, Vj) = 0. Para isto, 
observamos que 


A; (Vi, Vj) = pa 


; a 
/A)vj, pois À é simétrica 


| 
> 
as 
Pa 
Po 
« 
Re 


Portanto, (A; — À;) (vi,v;j) = 0. Como À;— À; £ 0, segue que 
(vi, Vj) = 0, isto é, os vetores v; e v; são ortogonais. [] 


O tipo de diagonalização que aparece no Exemplo 22.2 é muito im- 
portante na teoria das matrizes simétricas. Por isso, temos a seguinte 
definição. 
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Definição 22.2. 


Uma matriz A E M,(R) é dita diagonalizável por matriz ortogonal 
se existe uma matriz ortogonal P (lembre, P-! = Pº) e uma matriz 
diagonal D tais que A = PDP”. 


Da discussão do Exemplo 22.2, vimos que, para diagonalizar uma 
matriz A E M, (RR) utilizando uma matriz ortogonal P, foi preciso encon- 
trar n autovetores linearmente independentes e ortogonais. A questão 
é: quando é que isso é possível de ser realizado? O próximo teorema 
caracteriza o tipo de matriz que pode ser diagonalizada por matriz or- 
togonal. 


Teorema 22.3. 


Uma matriz 4 € M,(R) é diagonalizável por matriz ortogonal se e 
somente se 4 é uma matriz simétrica. 


Demonstração 


Uma das direções é muito simples de ser feita. Suponha que A seja 
diagonalizável por matriz ortogonal, como na Definição 22.2, então 


A! = (PDP) = (PY'D'P' = PDP'=A, 


onde (P)'=Pe D' =D, já que D é uma matriz diagonal. Assim, 
concluímos que 4 é uma matriz simétrica. 


A recíproca é muito mais complicada e será omitida nestas notas. A 
ideia básica desta parte da demonstração será apresentada na próxima 
aula e envolve um dos teoremas mais importantes da Algebra Linear. [] 


Exemplo 22.3. 


Determine se a matriz 


3 =. 4 
A=:| =20 6.52 
E 


é diagonalizável por matriz ortogonal e, caso seja, determine uma ma- 
triz ortogonal P e uma matriz diagonal D tal que A = PDP”. 


Solução: 


Como A é uma matriz simétrica, então, pelo Teorema 22.3, ela é diago- 
nalizável por matriz ortogonal. Vamos, agora, realizar o cálculo de diagonali- 
zação de 4. 


Os autovalores da matriz A são as raízes do polinômio característico 


p(x) = det(x—A) 


x—3 2 —4 
= 2 x—6 —2 
4 2 x—3 


= 02-12X2+421x+98. 


Observando, por inspeção, que À, = —2 é uma raiz de p(x), temos que 


plo = (x) (0 = 144449) = (44 D (x 7): 


Assim, os autovalores da matriz A são À, = —2, com multiplicidade algébrica 
1,e A = 7, com multiplicidade algébrica 2. 


Para o autovalor À, = —2, temos que os autovetores associados, v = (x,7,), 
satisfazem o sistema linear 


(—25 — A)v =0. 


Completando os cálculos, temos que o autoespaço E(—2) é um subespaço 
de dimensão 1 e é gerado pelo vetor v; = (—2,—1,2). 


Para o autovalor À, = 7, como já sabemos que a matriz A é diagonalizável, 
o autoespaço E (7) tem dimensão igual a 2. O fato interessante é que podemos 
construir uma base ortogonal de autovetores para esse subespaço E(7). Os 
autovetores v = (x, y, Z) associados ao autovalor À, = 7 satisfazem o sistema 
linear 
(7h -A)v=0. 


Usando as técnicas usuais para a resolução de sistemas lineares, obtemos 
que: 
E(7) = eR |Av=N) 
(veRê| (7h-A)v=0) 
= (xyz) CR | 2x+y-22=0). 


Para obter uma base ortogonal de E(7), observamos facilmente que 
vo = (1,0,1) E E(7). O outro vetor vs = (a,b,c) E E(7) deve satisfazer 
2a+b— 2c = 0 e ainda ser ortogonal a v2, isto é, (v2,V3) = 0, ou seja, a+c = 0. 
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Portanto, v3 = (a,b,c) deve satisfazer o sistema linear 


2a+b-2c=0 
a+c=0. 


Completando os cálculos, obtemos, por exemplo, vz = (—1,4,1). Observe 
que, pelo Teorema 22.2, o autovetor v; é ortogonal aos autovetores v> e v3, já 
que eles correspondem a autovalores distintos da matriz simétrica A. Assim, 
(vi, v2, V3 | é um conjunto ortogonal de autovetores da matriz A. Normalizando 
esses vetores, obtemos: 


rs SiS) 
Iva 
w=-2. =(1/V2,0,1/V7); 
Iva 
u= Til =(-1/VTB,4/VT8,1/VT8). 


Portanto, fu,,u>,u3) é uma base ortonormal de autovetores de 4. Com 
esses autovetores, obtemos a matriz P e com os autovalores, obtemos a 
matriz D: 


-2/3 1/2 —-1//18 200 
P=| -1/3 0 4//18 |; D=|[ 070), 
2/3 2/N2 1/v18 Quad 


de modo que A = PDP". 
Exercício 22.1. 


1. Mostre que se 4 é uma matriz simétrica, então A? também é uma 
matriz simétrica. 


2. Mostre que se A é uma matriz diagonalizável por matriz ortogonal 
então A? também é. 


3. Determine uma matriz ortogonal P e uma matriz diagonal D tal 
que A = PDP”, onde a matriz A é dada por 


|. do Ao sd) 
2 dW-d 
e Qu ml=2 
O 0-=2. 1 


Aula 2 3 


O TEOREMA ESPECTRAL 


Ao final desta aula, você deverá ser capaz de: 


compreender o significado do Teorema Espectral; 
compreender a decomposição espectral de matrizes si- 
métricas; 

aplicar os conceitos apresentados em exemplos impor- 
tantes. 


Álgebra Linear II | O Teorema Espectral 


Pré-requisitos 
Aulas 5 e 22. 
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O TEOREMA ESPECTRAL 


Nesta aula, continuaremos estudando as matrizes simétricas e fare- 
mos uma breve discussão do chamado Teorema Espectral para Matrizes 
Simétricas, mencionado na demonstração do Teorema 22.3 da aula pas- 
sada. Os detalhes da demonstração desse importante teorema serão 
omitidos nestas notas. Uma versão simples do Teorema Espectral é 
apresentada a seguir. 


Teorema 23.1 (Teorema Espectral para Matrizes Simétricas). 


Seja 4 E M,(R) uma matriz simétrica (isto é, 4º = A). Então vale: 


1. A matriz À possui n autovalores reais, contando suas multiplici- 
dades. 


2. A dimensão do autoespaço associado a cada autovalor À é igual à 
multiplicidade de À como raiz do polinômio característico de 4, 
isto é, a multiplicidade geométrica de À é igual à sua multiplici- 
dade algébrica. 


3. Os autoespaços são ortogonais entre si, isto é, os autovetores as- 
sociados a autovalores distintos são ortogonais. 


4. A matriz À é diagonalizável por matriz ortogonal, isto é, exis- 
tem uma matriz ortogonal P e uma matriz diagonal D tal que 
A = PDP. 


£3 Como já foi observado anteriormente, o polinômio característico 
de uma matriz A não possui necessariamente apenas raízes reais. 
Por exemplo, dada a matriz 


ds 0 —1 

(10) 

seu polinômio característico, dado por p(x) = x? + 1, não pos- 
sui raízes reais. Mas isso não acontece se A for uma matriz 
simétrica. O item 1 do Teorema Espectral afirma que o polinômio 
característico de uma matriz simétrica possui apenas raízes reais. 
A demonstração desse fato, embora simples, é bem trabalhosa 
e utiliza o Teorema Fundamental da Álgebra, que diz que todo 
polinômio de grau n com coeficientes reais possui n raízes reais 
ou complexas, contando suas multiplicidades. Na demonstração 


do Teorema Espectral mostra-se que as n raízes do polinômio ca- 
racterístico são, de fato, raízes reais. 

É Se A é uma matriz simétrica e tem n autovalores distintos, então 
pelo Teorema 5.2 da Aula 5 e pelo Teorema 22.2 da Aula 22, 
vemos que A é diagonalizável por matriz ortogonal. 


£ Se A é uma matriz simétrica e tem algum autovalor com multipli- 
cidade algébrica maior que 1, ainda é verdade que podemos di- 
agonalizá-la. Na verdade, podemos mostrar que se A é simétrica 
e tem um autovalor À de multiplicidade k, então o autoespaço 
associado tem dimensão k. Isto significa que o sistema linear 
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(Ah-—-A)v=0 


admite k soluções linearmente independentes, isto é, a matriz A 
tem k autovetores linearmente independentes associados ao auto- 
valor À. Usando o processo de ortogonalização de Gram-Schmidt, 
podemos obter uma base ortonormal para este autoespaço. Obte- 
mos assim um conjunto de k autovetores ortonormais associados 
ao autovalor À. Como autovetores associados a autovalores dis- 
tintos são ortogonais, então, considerando o conjunto de todos os 
autovalores de A, obtemos uma base ortonormal de autovetores 
para R?. Consequentemente, A é uma matriz diagonalizável, e a 
matriz diagonalizadora P, formada pela base de autovetores de A, 
é uma matriz ortogonal. 


DECOMPOSIÇÃO ESPECTRAL DE UMA MATRIZ 


SIMÉTRICA 
Seja A E M,(R) uma matriz simétrica e (u,,u2,...,u,) uma base 
ortonormal de autovetores associados aos autovalores Ay, Ã»,..., Ay da 


matriz 4. Seja P a matriz ortogonal tendo esses autovetores como colu- 
nas e D a matriz diagonal tal que A = PDP". Então 


A = PDP 
Ag dr age 20 
O: Ago se O 
= [lu uz ui) |. o uu cu 
O gs sas ado 
— [Mui Aous dad Al] [ui US e: u,] 


Muçu + Aquous + ---Auul,. 
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Esta representação é chamada uma decomposição espectral de A. 
Exemplo 23.1. 


Obtenha uma decomposição espectral da matriz A = ( , ; ) ; 


Solução: 


Sendo A uma matriz simétrica, essa decomposição existe. O polinômio 
característico de A é dado por 


p(x) = det(xb — A) 
= x-—1Ix+24 
= (x—-8)(x—3). 


Então os autovalores são Ay = 8 e À = 3, e ainda podemos obter os respec- 
tivos autovetores uy = (2/45,1/v5) eu, = (-1/v5,2/v'5). Assim, temos 
que 


A = PDP 
2 — 2 1 
72) | v5 v5|/80)[ 5 v5 
a) 12 Es E RR. 
v5 5 v5 5 


Denotando a matriz P = [u; us], temos, pela decomposição espectral, que: 


A =8uuí +3uçuí. 


Para verificar essa decomposição da matriz 4, observe que: 


2 dd 
v5 2 a o 
Rc RE (5% =) - dl 
is ss 

= e 52 

wu, = v5 (= =) - o 

2 v5 5 SD A 

v5 5 5 


e, finalmente, 


2a (3 6) (72 
5 5 5 
8uju/ + 3uu) — + = =A 
ria 6 8 6 12 2 4 
5 5 5 5 


PROCESSO DE DIAGONALIZAÇÃO DE UMA MATRIZ 
SIMÉTRICA 4 E M,(R) 


1º Passo: Obtenha o polinômio característico da matriz 4, 


p(x) = det(xl, — A). 


2º Passo: Encontre as raízes do polinômio característico de 4. Elas 
são todas reais e existem exatamente n delas, contando suas multiplici- 
dades. 


3º Passo: Para cada autovalor À da matriz A, de multiplicidade 
algébrica k, determine seu autoespaço associado 


E(A)=(veR"|(AL—A)v=0h, 


que é um subespaço vetorial de dimensão k. Para cada E(A) assim 
obtido, determine uma base ortonormal que consistirá de k autovetores. 


Se desejar, pode utilizar o processo de ortogonalização de Gram-Schmidt. 


A reunião dessas bases determina uma base ortonormal de autovetores 
para Rº. 


4º Passo: Seja P a matriz cujas colunas são os n autovetores da 
base ortonormal de R” obtida no terceiro passo. Portanto, P é uma 
matriz ortogonal. Seja D a matriz diagonal cuja diagonal principal é 
formada pelos n autovalores da matriz 4, tomados na mesma ordem de 
seus autovetores correspondentes na matriz P. Temos, então, 


A = PDP". 


Exemplo 23.2. 


Aplique o processo de diagonalização acima à matriz 
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e obtenha sua decomposição espectral. 
Solução: 


Observe, inicialmente, que 4 é uma matriz simétrica e, portanto, se aplica 
o processo de diagonalização acima. Não é difícil determinar que o polinômio 
característico da matriz A é dado por 


p(x) = det(xl; —A) = (x+2)*(x— 4), 
de modo que os autovalores de A são: 
AM =—2. com multiplicidade algébrica 2, e 


Ap =4 com multiplicidade algébrica 1. 


O autoespaço associado a À = —2 é dado por 


E(-2) (ve Rê | (A+26b)v=0) 


= ((xyp)eR|x+y+2=0). 


Para escolhermos uma base ortogonal de E(—2), podemos usar o pro- 
cesso de ortogonalização de Gram-Schmidt a partir de uma base qualquer 
de E(—2) ou podemos tentar obter diretamente dois vetores ortonormais de 
E(—2), como já foi feito anteriormente. Faremos o cálculo diretamente. Da 
equação x + y + z = O podemos ver facilmente que v; = (1,0,-1) E E(—2). 
O outro vetor, v> = (a,b,c) E E(—2), deve satisfazer a +b +c = 0 e ainda ser 
ortogonal a vi, isto é, (v2,v1) = 0, ou seja, a— c = 0. Portanto, v> = (a,b,c) 
deve satisfazer o sistema linear 


a+b+c=0 
a-c=0. 


Completando os cálculos, obtemos, por exemplo, vz = (1, 2,1). Norma- 
lizando esses dois vetores, obtemos: 


u = il + (1/V2,0,-1//D) e 
1 
m= = (1/V8,-2/Vô, 1/5). 


Assim, (uy, u>+ forma uma base ortonormal do autoespaço E (—2). 
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Por outro lado, o autoespaço associado a À, = 4 é dado por 


E(4) (veRº*| (46 -A)v=0) 


(yz) cR'|x=ze y=2). 


É fácil ver que v3 = (1,1,1) € E(4). Normalizando esse vetor, obtemos 
que 


u = Tl = (1/03,1/V3,1/43) 
representa uma base ortonormal do autoespaço E (4). Como A é matriz simé- 
trica, os autovetores associados a autovalores distintos são ortogonais e, assim, 
us é ortogonal a uj e uz. Portanto, (u,,u2,u3) é uma base ortonormal de Rº 
formada por autovetores de 4. Com esses autovetores, obtemos a matriz P e, 
com os autovalores obtemos a matriz D: 


1/2. 1/46 1/3 
P=[u, wu us])= 0 -2//6 1/3 |: 
1/02 1/06 1/V3 


EA TOO 
D=| 0 20, 
O 04 


de modo que A = PDP". A decomposição espectral da matriz A é dada por: 


A = —2uju — 2uzus + S4usus, 
ou ainda, e tm | 
3 0 —3 S 6. ao a 
Re RR E nO GPRS 
“gu O q E E. 6 Ed 
e gd a dd 
im cipa Goto aa 
E) e a 
022 
a ta 
220 
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Resumo 


É muito importante que você entenda bem o significado deste Teo- 
rema Espectral. Lembre do que aconteceu em exemplos vistos an- 
teriormente, em que a matriz considerada não era simétrica. Estu- 
damos exemplos de matrizes não-simétricas com autovalores repeti- 
dos que eram diagonalizáveis e outros exemplos de matrizes não- 
simétricas que não eram diagonalizáveis. Há algumas diferenças 
marcantes entre os casos simétrico e não-simétrico que tentaremos 
resumir agora. 

Se A for uma matriz não-simétrica, então nem todas as raízes de 
seu polinômio característico precisam ser números reais, o que é 
necessário no caso de a matriz A ser simétrica. Se A for uma matriz 
não-simétrica e todas as raízes de seu polinômio característico forem 
números reais, então ainda é possível que 4 não seja diagonalizável. 
É o caso em que um autovalor À de multiplicidade algébrica k não 
possui k autovetores linearmente independentes, isto é, quando o 
autoespaço correspondente tem dimensão menor que k, ou ainda, 
quando a multiplicidade geométrica do autovalor é menor que sua 
multiplicidade algébrica. Agora, quando A é uma matriz simétrica, 
além de todos os autovalores serem reais, são iguais a multiplicidade 
algébrica e a multiplicidade geométrica de cada autovalor. 

E, por fim, diferente do que ocorre no caso de matriz simétrica, se a 
matriz A é não-simétrica, então autovetores associados a autovalores 
distintos não precisam ser ortogonais. Estude e analise, com a ajuda 
de seu tutor, exemplos já vistos em aulas anteriores em que ocorrem 
as diferenças descritas aqui. 


Exercício 23.1. 


1. Em cada caso, aplique o processo de diagonalização à matriz 4, 
determinando matrizes ortogonal P e diagonal D, tais que 


A = PDP". 
2 
a A= ( 29 ) 
0 —1 —1 
bA4A=| 1 0 —1 
-1 -1 0 
2 2 0.0 
2 2.00 
CR QUO: rd 
00272 
é a RR | —1 
2. SejjamÃ=| 13 1 |ev= 1 |. Verifique que À = 5 
113 0 


é um autovalor de A e que v é um autovetor de 4. Em seguida 
obtenha matrizes ortogonal P e diagonal D, tais que A = PDP". 
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Aula D4 : 


OPERADORES AUTOADJUNTOS 


Ao final desta aula, você deverá ser capaz de: 


compreender o conceito de operador autoadjunto; 


aplicar os conceitos apresentados em exemplos impor- 
tantes. 


Álgebra Linear II | Operadores Autoadjuntos 


Pré-requisitos 
Aulas 8e 20a 23. 
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OPERADORES AUTOADJUNTOS 


Nesta aula vamos definir os operadores lineares T : IR” — R” asso- 
ciados às matrizes simétricas e estudar suas propriedades. Como estare- 
mos trabalhando sempre com bases ortonormais, é de suma importância 
que o espaço vetorial R” esteja munido de um produto interno, o qual 
estaremos sempre supondo que seja o produto interno canônico de R”. 


Definição 24.1. 


Um operador linear T : R”º — R” é denominado autoadjunto se sa- 
tisfaz 


(T(u),v) = (u,T(v)) para todo u,v e R”. 


O resultado que segue relaciona os operadores autoadjuntos com as 
matrizes simétricas. 


Teorema 24.1. 


Um operador linear T : Rº — R” é autoadjunto se e somente se a 
matriz 4, que representa T com respeito a qualquer base ortonormal a 
de R”, é uma matriz simétrica. 


Demonstração 


Com respeito à base ortonormal a de R”, temos que T (u) = Au para 
todo u E R”. Assim, para todo u, v E R”, temos que 


(Tu,v) = (Au,v) = (Au)'v = u/A'y 


(u, Tv) = (u,Av) = u'Ay, 


onde A! é a transposta da matriz A. Assim, 


T é autoadjunto <> (T(u),v) = (u,T(v)) para todo u, v ER? 
(Au, v) = (u,Av) para todo u,v E R? 
u'A'y = u'Av paratodou, ve R” 

AÍ=A 

A é uma matriz simétrica . 


t1%4 


E importante salientar que não existe uma relação tão simples entre 


o operador linear T : R” — R” e sua representação matricial A = [T]w 
quando a base o não for ortonormal (veja a observação ao final do 
Exemplo 24.1). 


O Teorema 24.1 também fornece um critério prático para determinar 
se um dado operador linear T : R” — R” é autoadjunto. Basta considerar 
qualquer base ortonormal à de R” e verificar se a matriz A = [T|o é uma 
matriz simétrica. 


Exemplo 24.1. 


Determine se o operador linear 


é autoadjunto. 
Solução: 


Vimos, no Exemplo 20.1 da Aula 20, que T é a projeção ortogonal sobre o 
eixo-x. Considerando a base canônica q = fe,,e>+ de R2, vimos que a matriz 
que representa T' nesta base é dada por 


a=Iro=(0 oi 


Como a base canônica é ortonormal e a matriz A é simétrica, então, pelo 
Teorema 24.1, o operador T é autoadjunto. 


Vejamos o que acontece quando escolhemos um base 8 de R? que não é 
ortonormal. Considere a base 8 = (u,,u>+ dada por 


u = (V2/2,V2/2) e w = (0,1). 


Está claro que esta base não é ortonormal, e ainda temos que 


Tu =T(V3/2,3/2) = (42/2,0) = 1-m + (=3/2)-uo 
Tu, =7(0,1) — (0,0) =0-u,+4+0-u,. 


Daí, segue que a matriz que representa T na base f é dada por 


s=m=( 5h 0) 


Observe que esta matriz não é simétrica, mas também a base 2 não é ortonor- 
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mal, o que não contradiz o Teorema 24.1. 


Exemplo 24.2. 


Considere os operadores lineares 


T:RESR?, Ti(x,)) = (2,27) 


Bb: Rº — Rê D(x,y) = (3,2). 


Verifique que 74 e 75 são operadores autoadjuntos e verifique se a 
composição 71 o T) também é operador autoadjunto. 


Solução: 


Considerando a base canônica de R?, verificamos que as matrizes A, e A 
que representam respectivamente, os operadores Tj e T> nesta base, são dadas 


por 
1 0 01 
ds pes 
Como essas duas matrizes são matrizes simétricas, concluímos, pelo Teo- 


rema 24.1, que T1 e T» são operadores autoadjuntos. No entanto, o operador 
obtido pela composição 


Toh: R? ==> Rº, (n 0 T)(x,)) = (7,2x) 


é representado, na base canônica, pela matriz 


01 
»=(2 0): 


que não é uma matriz simétrica. Assim, outra vez pelo Teorema 24.1, a 
composição T| o T não é um operador autoadjunto. Daí, concluímos que a 
composição de operadores autoadjuntos não é, necessariamente, autoadjunto. 


O próximo teorema segue imediatamente dos resultados sobre ma- 
trizes simétricas estudados nas Aulas 22 e 23. 


Teorema 24.2. 


Seja T:R” — R” um operador autoadjunto. Então 


1. Autovetores correspondentes a autovalores distintos de T são or- 


togonais, isto é, se V1,V2,..., Vk São k autovetores associados aos 
autovalores distintos A, Ã»,..., Ay, então vi,V>,...,Yy São orto- 
gonais. 


2. O operador T possui n autovalores reais, contando suas multipli- 
cidades. 


3. A dimensão do autoespaço associado a cada autovalor À é igual à 
multiplicidade de À como raiz do polinômio característico de T, 
isto é, a multiplicidade geométrica de cada autovalor À é igual à 
sua multiplicidade algébrica. 


4. Os autoespaços de T são ortogonais entre si. 


5. Existe uma base ortonormal (u,,us,...,u,! de R” formada por 
autovetores de 7. 


A última afirmação do Teorema 24.2 também é conhecida como 
Teorema Espectral para Operadores Auto Adjuntos Reais e diz, simples- 
mente, que estes operadores são diagonalizáveis. 


Exemplo 24.3. 


Seja T : Rê — IR? dado por 


T(x,7,2) = (3x,2y +2,y +22). 


a. Verifique que T é um operador autoadjunto. 
b. Determine os autovalores e os autovetores de T e verifique que T 
é diagonalizável. 


Solução: 


a. Considerando a base canônica fe,,e>,e3) de Rº, temos que 


Te; =T(1,0,0) = (3,0,0), 
Teo =T(0,1,0) = (0,2,1), 
Tes=T(0,0,1) = (0,1,2). 
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Assim, a matriz que representa o operador linear T na base canônica é 
dada por 


3 00 
A=| 021 
012 


Observando que A é uma matriz simétrica, temos, pelo Teorema 24.1, que T é 
um operador autoadjunto. 


b. O polinômio característico do operador T é dado por 


p(x) = det(xl;— A) 
x—-3 0 0 
= O x—-2 + 
0 =| 2 
= (x-3)2(x—1) 


Assim, os autovalores de T são A = 3, com multiplicidade algébrica 2, e 
À» = 1 com multiplicidade algébrica 1. Não é difícil obter que o autoespaço 
E(3), associado a À = 3, é dado por 


E(3) = fveR|Tv=3v) 
= [((x,y,2) €R'|y=ze x arbitrário). 


Portanto, uma base ortonormal de E (3) é dada por 


Std ss (0.55) 


Analogamente, o autoespaço E (1), associado a À> = 1, é dado por 


E(1) = fveR| Tv=v) 
(xp) eR|x=0ey=-z), 


1 —1 
e uma base ortonormal de E(1) é dada pelo vetor us = (o. NG 3) Con- 
sequentemente, 8 — fu,,u>,u3! é uma base ortonormal de IR? formada por 
autovetores de T e, nesta base, T é representado pela matriz diagonal 


300 
B=[7M]p=[ 030 
001 


Portanto, T é um operador diagonalizável. 
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Exemplo 24.4. 


Determine valores de a, b € R de modo que o operador T:Rº — Rº, 
definido por 


T(x,7,2) = (x+2ay +27,4x — 5y — bz, 2x — 4y +27), 


seja autoadjunto. Determine, também, uma base ortonormal de Rê for- 
mada por autovetores de T e a matriz que representa T nesta base. 
Solução: 


Considerando a base canônica [e;,e>,e3) de IRÊ, temos que 


Te =T(1,0,0)=(1,4,2)=1.e;+4-e+2-e5, 
Te,=T(0,1,0)=(2a,-5,-4) =2a-e;+(—5)-e +(—4) -es, 
Tes = T(0,0,1) = (2,—b,1) =2-e,+ (—b) "e+l-es. 


Assim, a matriz que representa o operador linear T na base canônica é 
dada por 


1 2a 2 
A=| 4 -5 —b 
2 —4 1 


Para que T seja um operador autoadjunto, é necessário que a matriz A seja 
simétrica, isto é, que 4 = A. Para isso, é preciso que 2a =4e —b= —4, ou 
seja, que 

a=2eb=4. 


Assim, obtemos a matriz simétrica 


| ÃO 
A=| 4 -5 4 |, 
7 | 


garantindo que o operador T é autoadjunto. Não é difícil verificar que o 
polinômio característico de T é dado por 


det(x15 —A) 
(x +9)(x — 3)2. 


p(x) 


Os autoespaços correspondentes são dados por 


E(-9) = (veRê| Tv=-9) 
= [(gyp)ceR|x=-zey=2), 
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E(3) = (veRº| Tv=3v) 
= (no) eR'| -x+2+7=0). 
Uma base ortonormal de E (—9) é dada pelo vetor u, = (ve Te)» en- 
quanto uma base ortonormal de E (3) é dada pelos vetores uz = (5 0, 5) 


DO do EN 1 Es ; 
e us = (5 cl a)» Consequentemente, B = (u,,u>,u3) é uma base 


ortonormal de R? formada por autovetores de T e, nessa base ordenada, T 
é representado pela matriz diagonal 


-9 0 0 
B=[T], = 030 
003 
Observe que T é um operador diagonalizável. 


Exemplo 24.5. 


Dados os vetoresu = (4,4, -2),v=(4, 2, 4)ew=(1, —2, —2), 
seja T : IRº — Rº o operador linear dado por 


Tu=(10,-2,-2), Tv=(-2,10,-2) e Tw= (1,1,-5). 


Verifique que T é um operador autoadjunto. 
Solução: 


E fácil ver que (u, v, w+ é uma base ortogonal, pois 


(u,v) =4:444.(-2)4(-2).4=0; 
(u,w) =4-:1+4:(-D+(-2)-(-2)=0; 
(vw) =4:1+(-9-(-D)+4-(-2)=0 


Assim, os vetores normalizados 


u= fal = na, 
w= Te (2/3,-1/3,2/3) e 
w 
w==(1/3,-2/3,-2/3) 
ud] 
formam uma base ortonormal de Rê. Como |[ul| = ||v|| = 6 e ||w]|| = 3, temos 
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T(u,) a) x r(gu) E ET (u) = 1(10,-2,-2)= 
(1553): 
T(u,) (o) (5 ) E ET(v) E 5(-2,10,-2) E 
—1l 5 —1 
- (555): 
T(us) = (o) =r(2w) = ST() = (LL-5)= 


Agora, não é difícil ver que os vetores T(u,),T (u>) e T (us) se expressam 
em função da base 8 = (u,,u>,u3+ como: 


T(u;))=(5/3,-1/3,-1/3)=1-u+1-w+1-us; 
T(u,s)=(—1/3,5/3,-1/3) = 1:u + (1) -u + (—1)-us; 
T(us) = (1/3,1/3,-5/3)=1-u +(—1)-u,+1-us. 


Portanto, a matriz que representa o operador T com respeito à base ortonor- 
mal (u,,u>,us) é dada por 


Ee dó À 
B=[Tp=[ 1-1 =1 
UR | 


Como B é uma matriz simétrica, concluímos, pelo Teorema 24.1, que o 
operador T é autoadjunto. Observe que, neste exemplo, usamos uma base 
ortonormal que não é a base canônica nem é uma base de autovetores. 


Autoavaliação 


E de suma importância que você reveja e entenda muito bem a 
relação que existe entre as matrizes simétricas, estudadas nas aulas 


anteriores, e os operadores autoadjuntos vistos nesta aula. Compare 
os conceitos e estude os exemplos. Em caso de dúvidas, não hesite 
em consultar o seu tutor. 
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Exercício 24.1. 


1. Verifique que o operador T : 


Rê 


Rê, dado por 


Ty) — (2x+y+z, Pay E: x—y+22), 


é autoadjunto. 


2. Determine uma base ortonormal de autovetores do operador T 


dado no exercício anterior. 
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Aula 2 b 


FORMAS BILINEARES 


Ao final desta aula, você deverá ser capaz de: 


compreender o conceito de forma bilinear; 


aplicar os conceitos apresentados em casos particula- 
res. 
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FORMAS BILINEARES 

Pré-requisito 

Aula 22. Nesta aula vamos introduzir um conceito que generaliza a noção 
de aplicação linear num espaço vetorial. Mais especificamente, va- 
mos desenvolver o conceito de forma bilinear, que dá origem às formas 
quadráticas que serão estudadas na próxima aula. Veremos a definição 
de formas bilineares e estudaremos algumas de suas propriedades, prin- 
cipalmente sua relação com as matrizes, o que constitui o aspecto mais 
importante para fins práticos. 


Definição 25.1. 


Seja V um espaço vetorial real. Uma forma bilinear em V é uma 
aplicação 
B:iVxV 5 R 
(u,v) » B(u,v) 


que é linear em cada uma das duas variáveis u e v, isto é, que satis- 
faz: 


1. paratodou, v weVeaeR, 


B(u+w, v)=B(u, v) +B(w, v) 
B(au, v)=aB(u, v); 


ii. paratodou, v weVeaceR, 


B(u, w+v)=B(u, w+B(u, v) 
B(u,av)=aB(u, v). 


Exemplo 25.1. 


Seja F o produto escalar em V = Rº, isto é, dados 
u=(m,u2,.. Un), V=(vi,v2,..., Vn) ERº, considere a aplicação 


F:VxV Ss R 
(u,v) +» F(u,v)=uvi+uvo + und. 


Verifique que F é uma forma bilinear em R”. 
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Solução: 

De fato, considerando outro vetor w= (Wi, W2,..., vn ER" caeR, 
temos que 
F(utaw,v) = B((wm-+Haw, uz taw,..., Untawn), (Vi, V2,..., Vn)) 


(um +awy)vi + (uz +awo)vo +---+ (Un +awn)Vn 
— (uvituvo + +unva) +Ha(wivi + wav +: + Wav) 
= F(u,v)+aF(w,yv), 


o que mostra que F(u, v) é uma transformação linear na primeira variável 
u. Um argumento análogo, deixado a cargo do aluno, mostra que F(u, v) 
também é uma transformação linear na segunda variável v. Assim, podemos 
concluir que F(u, v) é uma aplicação bilinear de R”. 


Exemplo 25.2. 


Seja a matriz 
200 
A=|420 
003 


N 


Mostre que podemos associar à matriz A uma forma bilinear 
B:RºxRº — R dada por 


2 Ag MO » 

B((x1, x2,x3), (01, 72,73) = (01x02x3)| 420 X 
003 y. 

= 2x1714+42x271 +2x272 431373 


Solução: 


Observe que para todo par de vetores u, v c Rº 


X1 y1 
u= X2 eVv=— 2 |» 
X3 3 


podemos reescrever 
B(u, v) =u'Av, 


onde u” é a matriz transposta de u. Assim, a bilinearidade da aplicação B(u, v) 
decorre facilmente das propriedades do produto e da soma de matrizes. 


Este exemplo é facilmente generalizado. 
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Teorema 25.1. 


Seja A = (a;;) E M,(R), isto é, uma matriz de ordem n. Podemos 
associar à matriz A uma forma bilinear F :R” x IR” —» R dada por 


F(u, v) =u'Ayv, 


ondeu, ve R”. 


Observe que, reescrevendo os vetores u e v na forma 


X1 X1 
X2 y2 
u= ev= ; , 
Xn Yn 
então 
F(u,v) =u'Av 
dj dj co din X, 
dj 42 e dIn » 
=(X1 X2 *** Xa) A 
nl An tt Am A 
= 011411 + 012X1)2 + *** + AnnXnYn 
n 
=) aja; 
ij=1 


Sejam V um espaço vetorial real, F:V xV — R uma forma bilinear em 
V,ea = (e;,e>,..., e! umabase de V. Sejamu, ve V com 


u=ujej+ue+:--+ une 


V=vVi€ +ve +: +Vnen. 
Então, 
F(u,v) =F(ue+ue+---+u,en, Vier +vwe +: + van) 
= ujviF(e,,e1) + ujv2F (e1,e2) Soa +unVnF (en, En) 
n 
= D uiv;F(e;,€;) ê 


ij= 


Assim, a forma bilinear F fica completamente determinada pelos n? 
valores F(e;, e ;). 


Definição 25.2. 


A matriz A = (a;;), com a;; = F(e;, e;), é chamada de representação 
matricial da forma bilinear F com relação à base a, ou, simples- 
mente, de matriz de F com relação a o. 
Esta matriz representa F no sentido que 


n 
= umjie e = 


Pad 


= [ulo4 [va 


para todo par de vetores u, v e V. Como de costume, [ul]; denota o 


vetor das coordenadas de u com respeito à base o. 


Exemplo 25.3. 


Seja a forma bilinear F : 


RZ xR? > 


R dada por 


F(u, v)= F((x1, x2), (y1, )2)) = x171 — 2172 4320271 — 5 x2y2, 


para todo u = (x1, x2), v = (y1, y2) E R2. Considere a = fey, esj a 
base canônica de R? e 8 = ((1, 0), (1, 1)J outra base de R?. Deter- 
mine a matriz de F com respeito a essas bases. 


Solução: 


Primeiramente, façamos o cálculo da matriz de F com respeito à base 
canônica: 


F(er, e1) =F((I, 0), (1, 0)) A 
F(e,, e2) =F((I, 0), (0, D) Ei O 
F(e>,em)=F((0, 1), (1,0) =3; 
F(es, e2) F((o, D, (0, D) usa 


Portanto, temos que a matriz de F na base canônica é 


ER 


Para a matriz de F na base B, temos 


PULO (140))=1; 
F((1,0), (1,1))=0; 
F((1,1), (1,0))=4; 
FL, A (1, 1))=-2. 


Portanto, temos que a matriz de F na base 8 = ((1,0), (1, 1)+é 
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1 0 
ge 
Um problema interessante é saber qual a relação entre as matrizes 4 e B que 


representam uma mesma forma bilinear F em duas bases a e B, respectiva- 
mente. 


No caso do exemplo anterior, se P representa a matriz mudança de base, 
da base B para a base a, temos 


Daí, 


| O 
ae 
=P'AP. 
De um modo geral, temos o seguinte teorema: 


Teorema 25.2. 


Seja F uma forma bilinear de um espaço vetorial V. Se À é a matriz 
de F numa base & e B é matriz de F numa base B de V, então 


B=P'AP, 


onde P é a matriz mudança de base, da base B para a base a. 


Definição 25.3. 


Uma forma bilinear F no espaço vetorial V é denominada simétrica 
se 
F(u,v)= F(v,u) 


para todo par de vetoresu, ve V. 


Teorema 25.3. 
Seja F uma forma bilinear no espaço vetorial V e 4 a matriz que 


representa F numa base a de V. Então F é uma forma bilinear simétrica 
se e somente se A é uma matriz simétrica. 
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Demonstração 


Por F ser uma forma bilinear em V, temos que 


= (u'Av), pois u'Av é um escalar 


Se, ainda, F for uma forma bilinear simétrica, então 
vAu=F(u,v)=F(v,u)=v'Au 
para todou, v E V. Portanto, temos 
AESA, 


isto é, a matriz A é simétrica. 


Reciprocamente, se A é uma matriz simétrica (isto é, 4 = A), então 


a forma bilinear F também é simétrica, pois 


F(u,v) =u'Av 
= (u'Av), pois u'Av é um escalar 
=vA'u 
=vAu, pois A/=A 
=F(v,u) 


para todo par de vetoresu, ve V. 


Autoavaliação 


Você deve ter compreendido que o conceito de forma bilinear é uma 
generalização do conceito de transformação linear já bastante estu- 
dado. E de extrema importância rever todos os conceitos e tentar 


resolver os exercícios propostos. Caso surjam dificuldades, consulte 
as notas de aula ou peça ajuda ao seu tutor. Os conceitos desta aula 
ainda serão bastante utilizados. Por isso, não deixe de fazer uma boa 


revisão de matrizes simétricas. 
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Exercício 25.1. 


1. Seja A E M,(R). Verifique que a aplicação F:RºxRº 5R, 
definida por F(u, v) = u/Av é uma forma bilinear. 


2. SejaF:RºxRº >, definida por F(u, v) = (u, v), o produto 
escalar em Rº. 


a. Determine a matriz A que representa a forma bilinear F com 
respeito à base canônica a C Rº. 


b. Determine a matriz B que representa a forma bilinear F com 
respeito à base 


BAU On t=E 010,21. 


3. Seja a forma bilinear F : R? x R? — R definida por 


F(u, v)=F((x1, x2), (71, 72) = 2x171 — 3 X172 +22, 


para todou = (x1,x2), v= (71, 2) ER2. 
a. Determine a matriz A que representa F com respeito à base 
ou =((1,0), (1, 1)+. 
b. Determine a matriz B que representa F com respeito à base 
B = (2, 1), Eine via 
c. Determine a matriz mudança de base P, da base B para a 
base q, e verifique que B = P'AP. 
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Aula 26 : 


FORMAS QUADRÁTICAS 


Ao final desta aula, você deverá ser capaz de: 


compreender o conceito de forma quadrática; 


aplicar os conceitos apresentados em casos particu- 
lares. 


Álgebra Linear II | Formas Quadráticas 


Pré-requisitos 
Aulas 22 e 25. 
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FORMAS QUADRÁTICAS 


As formas bilineares, vistas na aula anterior, dão origem às formas 
quadráticas que serão estudadas nesta aula. As formas quadráticas ocor- 
rem com grande destaque em aplicações da Álgebra Linear à Enge- 
nharia, como em critérios para projetos, em problemas de otimização 
e em processamento de sinais. Elas também ocorrem na Física, em 
descrições de energia potencial e energia cinética; em Economia, nas 
funções de utilidade; e, também, em Estatística. Em todas essas situ- 
ações é muito importante o conhecimento do sinal (positivo ou nega- 
tivo) que a forma quadrática pode assumir, assim como o conhecimento 
de seus autovalores associados. Uma parte muito importante da base 
matemática para o estudo das formas quadráticas segue facilmente do 
nosso estudo prévio sobre matrizes simétricas. 


Definição 26.1. 


Seja V um espaço vetorial real. Uma aplicação gq: V5R é 
chamada de forma quadrática se existe uma forma bilinear simétrica 
F:VxV>Rtalqueg(v) =F(v, v) para todo v ev. 


Seja A a matriz que representa a forma bilinear F na base à € V. 
Dizemos que matriz 4 é a representação matricial da forma quadrática 
q com respeito a essa mesma base q c V. Como a forma bilinear 
F é simétrica, então, pelo Teorema 25.3, da Aula 25, a matriz A é 
uma matriz simétrica. Com respeito à base o, denotamos A = (a;;) e 
v=(x1,x2;..., Xn) EV; então 


q(v) =F(v,v) 


= vAv 
am a tt An X1 
aj add An X2 
É X2 ... Xi) 
An1l And tt Am Xn 
n 
= D DijXiX j 
ij=1 


E agora, sendo A simétrica, vale que a;; = a ;i. Portanto, 


n n 
g(v) = D AijxiX; = anx Hans + Hama + 2) aijxix;. (26.1) 
ij i<j 


Observe ainda que, se A for uma matriz diagonal, isto é, a;; = O para 
i j, então teremos 


n 
Vaxix; — 0, 
i<j 
o que nos dá 
= 2 2 2; 
g(v) = aux +Hazgxs ++ anda, 


que será denominada representação diagonal da forma quadrática q. 
Veremos, mais à frente, que toda forma quadrática sempre admite uma 
representação diagonal. 


Exemplo 26.1. 


Seja a forma quadrática q : R? — R dada por 


g(x, y) = x — 10xy + 2. 


Determine a matriz A que representa a forma quadrática q com respeito 
à base canônica. 


Solução: 


Como A é uma matriz simétrica, podemos denotar 
a b 
RECO: 
a b bo 
q(x, )) (6 (5 2) E 


— ax? +2bxy + cy”. 


temos, então, 


Então, vale que 
a” +2bay+cy =x — 10xy +72, 


de onde concluímos que 


obtendo 
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(15) 


Observe que q é a forma quadrática associada à forma bilinear 


F(u,v) =(x 25) ( E E ) o ) 


=x1y1— 5%271 — 5x1)2 +X2)2, 


onde u = (x1, x2), Y= (y1, 2) E R?, com respeito à base canônica. 


Exemplo 26.2. 


Seja q: Rº — R a forma quadrática dada por 


g(v) = q(x1, X2, 1x3) = 5x7 +35 4225 — x1x2 +8x023, 


onde v = (x1, x2, x3) € Rº. Determinar a matriz A que representa a 
forma quadrática q com respeito à base canônica e expresse a forma 
quadrática na forma matricial g(v) = v'Av. 


Solução: 


Os coeficientes de x7, x3 e x formam a diagonal principal da matriz 4, 


como indica a equação (26.1). Como A é matriz simétrica, o coeficiente de 
X:X;, para i 4 j, é a soma dos coeficientes iguais a;; = a ji, como indica outra 
vez a equação (26.1). Portanto, 


1 
Mj=4p=5 (coeficiente de x;x;). 
Assim, é fácil ver que 
5 —-1/2 0 
A=| 1/2 3 4 
0 4.2 
E, finalmente, 
—1/2 0 X1 
g(x, x2,x3)=(x1 22 2x3) | —1I/2 0 3. 4 x2 
4 2 X3 


Queremos agora estudar o efeito de uma mudança de base sobre uma forma 
quadrática. Assim, sejam q: V — R uma forma quadrática e o e B duas bases 
do espaço vetorial V. Seja P a matriz mudança de base, da base B para a base 
a. Se À é a matriz que representa a forma quadrática q na base a e B é a matriz 


de q na base p, então, pelo Teorema 25.2, da Aula 25, sabemos que 
B=P'AP. 


Observe que, se P é uma matriz ortogonal que diagonaliza a matriz 4, então 
B=P'AP = P-lAP é uma matriz diagonal. Nesse caso, a matriz P também é 
chamada mudança de variáveis. Usaremos esses fatos no próximo exemplo. 


Exemplo 26.3. 


Determine uma mudança de variável P que transforma a forma qua- 
drática q: R? — RR, dada por 


2 é) 
g(x1, x2) = xj — 8x1x2 — 5x5 


na base canônica, em uma forma diagonal. Obtenha, também, a ex- 
pressão dessa forma diagonal. 


Solução: 


Observando os coeficientes de q, vemos que a matriz A que representa q 
na base canônica é dada por 


(455) 


Diagonalizar a forma quadrática q é equivalente a diagonalizar a matriz 
simétrica A. Usando os procedimentos já conhecidos sobre diagonalização 
de matrizes simétricas, os autovalores da matriz A são A —3e Ah, =-7. A 
matriz P será obtida a partir de uma base ortonormal de autovetores. Efetuando 
os cálculos, que é um exercício para você, obtemos 


w=( 2a 


, autovetor associado ao autovalor À =3,e 
-1/45 


w=( a 


, autovetor associado ao autovalor À, = —7. 
2/45 É 


Como (u,, u2) forma uma base ortonormal de R?, então 


INE ANE 
P=[u; ul=( is Rh 


e a matriz diagonal correspondente será 


30 
Be: 
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onde D = P'AP. 


A forma diagonal de q é dada por 


iss) 


onde 


é a mudança de variáveis. 
Veja que 

a(v) =q(x1,x2) =) — 8xyxp — 5x5 
E 1 —4 X1 
com ( 4 (5) 
=v'Av 
= (PwJ'A (Pw) 
=w'(P'AP)w 


=w'Dw 


Observe que a forma diagonal 
q(01, 2) =3yi— 1 
não contém o termo cruzado y1y2. 


Este exemplo anterior ilustra o teorema a seguir. A parte essencial 
de sua demonstração foi apresentada nos cálculos do Exemplo 26.3 e 
consiste na mudança de variáveis efetuada. 


Teorema 26.1 (Teorema dos Eixos Principais). 


Seja g:V — R uma forma quadrática. Então, sempre existe uma 
mudança de variáveis P que transforma a forma quadrática g(v) = v'Av 
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na forma diagonal g(w) = w' Dw, onde v = Pwe D=P'AP. 


O nome Teorema dos Eixos Principais segue do fato de que as co- 
lunas de P são chamadas eixos principais da forma quadrática q. Uma 
interpretação geométrica deste teorema será vista nas próximas aulas, 
mais precisamente no estudo da classificação de curvas cônicas e na 
classificação de superfícies quádricas. 


Exemplo 26.4. 


Determine uma mudança de variável P que transforme a forma qua- 
drática q: Rº — RR, dada por 


g(x1, X2, X3) = Es + 2a +25 + 4x1x2 + 4x02x3 


na base canônica, em uma forma diagonal. Obtenha também a ex- 
pressão dessa forma diagonal. 


Solução: 


Observando os coeficientes de q, vemos que a matriz A que representa q 
na base canônica é dada por 


> 

I| 
SIN QU 
[SS RS) 
= No 


Procedendo à diagonalização da matriz simétrica A, deixamos os detalhes 
dos cálculos como um exercício para você, obtemos os autovalores À, = 5, 


A =2eAy=-—l1. A matriz mudança de variável P será obtida a partir de uma 
base ortonormal de autovetores. Efetuando os cálculos, obtemos: 
2/3 
u,= | 2/3 | autovetor associado ao autovalor À, = 5; 
1/3 
—2/3 
Ww= 1/3 | autovetor associado ao autovalor À, = 2; 
2/3 
173 
u3= | —2/3 | autovetor associado ao autovalor Às = —1. 
2/3 
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Como (u,, u2, us) forma uma base ortonormal de Rº, então 


ja 9% AB 
P=[u; 5) us] = 213 1/3 —2/3 
1/3 2/3 2/3 


é uma matriz ortogonal e a matriz diagonal correspondente será 


5 0 0 
D=|/0 2 0], 
00 —1 
onde D = P'AP. 
A forma diagonal de q é dada por 
50 0 y1 
901,7») =017293)| 02 0 »2 
0 0 —1 3 
= 5yj +25 33, 
onde 
X1 1 
v= | x» ew=| 2 |; 
X3 Y3 
e 


v=Pw, ou w=P'y 
é a mudança de variáveis requerida. 
Observe, mais uma vez, que a forma diagonal 
q01, 32,3) = 51 +25 —)5 


não contém os termos cruzados y1y2, Y1Y3 € Y2Y3, isto é, os termos y;y; com 


iz j. 
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Aula ) 7 


CÔNICAS 


Ao final desta aula, você deverá ser capaz de: 


compreender o conceito de cônica; 


aplicar os conceitos apresentados em casos particu- 
lares. 


Álgebra Linear II | Cônicas 


Pré-requisitos 
Aulas 22, 25 e 26. 
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CÔNICAS 


Nesta aula estudaremos algumas figuras importantes do R?, ou seja, 
determinados conjuntos de pontos do plano cujas coordenadas satisfa- 
zem certas propriedades. Mais precisamente, consideraremos subcon- 
juntos de R? cujas coordenadas (x, y) satisfazem uma equação do tipo 


a? +by+cy) +dxrey+f=o, 


onde a, b, c, d, ee f são constantes reais (com pelo menos um dos 
números a, b ou c diferente de zero). A ideia toda é simplificar e clas- 
sificar equações desse tipo e, para isso, usaremos os resultados sobre 
diagonalização de formas quadráticas apresentados na aula anterior. 


Definição 27.1. 


Uma cônica é um conjunto de pontos do R? cujas coordenadas (x, y), 
em relação à base canônica, satisfazem uma equação do tipo 


ao +hy+ey +ixto+f=o, (27.1) 


onde os coeficientes a, b, c, d, e e f são números reais e pelo menos 
um dos números a, b ou c é não-nulo. 


Observe que a equação (27.1) contém uma forma quadrática, 
gx, )=a+by+oy”, 


uma forma linear, 
f(x,y) = dx+ey, 


e o termo constante f. 


Exemplo 27.1. 


Identifique o conjunto dos pontos (x, y) € R? que satisfazem a e- 
quação 
pre y —4=0. 


Solução: 
Comparando a equação 
do) 


com a equação (27.1), vemos que os valores dos coeficientes são a = c = 1, 
b=d=e=0ef= —4,e, portanto, representa uma cônica. Reescrevendo a 
equação na forma 


R+ry=A, 


identificamos os pontos (x, y) como pertencendo à circunferência de centro 
(O, 0) e raio 2, como ilustra a Figura 27.1. 


Figura 27.1: A circunferência x? +y? = 4. 


Exemplo 27.2. 


Identifique o conjunto dos pontos (x, y) € R? que satisfazem a e- 
quação 
y — kx=0, 


onde k é um número real não-nulo. 
Solução: 


Comparando a equação 
y-kx=0 


com a equação (27.1), vemos que os valores dos coeficientes são c = 1, 
a=b=e=f=0ed=-k+0,e, portanto, representa uma cônica. Re- 
escrevendo a equação na forma 


rem 
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identificamos os pontos (x, y) como pertencendo a uma parábola com eixo 
coincidindo com o eixo-y, como ilustra a Figura 27.2. 


YA X, 
0 X 0 bg 
caso k<Q caso k>0 


Figura 27.2: A parábola y? — kx. 


Exemplo 27.3. 


Identifique o conjunto dos pontos (x, y) € R? que satisfazem a e- 


quação 
x y e 
a bo o 

coma, beR,a, b>o0. 

Solução: 

Comparando a equação 

2 y 
pe 
o db 


com a equação (27.1), vemos que ela também representa uma cônica. Rees- 
crevendo a equação na forma 


a o x 
bB q? 
temos 
b 
À = E é 


o que representa um par de retas concorrentes que passa pela origem, como 
ilustra a Figura 27.3. 
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Figura 27.3: Asretasy = +2x. 


Os próximos exemplos mostram como procedemos para simplificar 
uma equação de uma cônica. 


Exemplo 27.4. 


Identifique a cônica representada pela equação 


5xº —4xy+8y) —36=0. 


Solução: 


Precisamos, inicialmente, eliminar o termo misto (—4xy); para isto, reali- 
zamos diagonalização da forma quadrática correspondente, 


g(x, y)=5xº —4xy + 8y2. 
Escrevemos a equação 5x? — 4xy + 8y? — 36 = 0 na forma matricial 


vAv=36, 


= X 2 o 5 —2 
(5 )er a e 


Lembre-se de que, na Aula 26, a matriz A é a matriz simétrica que re- 
presenta a forma quadrática g(x, y) = 5x? — 4xy + 8y? com respeito à base 
canônica. Não é difícil ver que os autovalores da matriz 4 são À; =4e A =9, 
e os autovetores normalizados são 


wa 


autovetor associado ao autovalor Ày = 4 
1 / VS ) ) 1 


com 
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w=( Rio 


autovetor associado ao autovalor À, = 9. 


Como (u;, u2) forma uma base ortonormal de R?, então 


ENC E o 
P=imml=( 11 2/8) 


é uma matriz ortogonal que diagonaliza a matriz A e a matriz diagonal corres- 


pondente será 
40 
Belo 


A forma diagonal de q é dada por 


atum =t) (0 9) (5) 


= 4x? + 9y% R 


(5) e nelã): 


v=Pv,, ou vi=P'y. 


Temos que D = P'AP. 


onde 


com 


Portanto, a equação da cônica pode ser reescrita como 


ou ainda, 
dx +9y7 = 36, 


o que nos dá a equação 
MM 
(li + Es 
9 4 
que representa uma elipse de semieixo maior 3 e semieixo menor 2, como 


ilustra a Figura 27.4. 


=], 
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2 2 
Figura 27.4: A elipse | + 1 =1. 


Exemplo 27.5. 


Identifique a cônica representada pela equação 


22 + 4x) +27) +4V2x+12/2y—-8=0. 


Solução: 


Observe que neste exemplo a forma linear /(x, y) = dx+ey = 
42x + 1242 y é não-nula. Reescrevendo a cônica na forma matricial, obte- 
mos 


vVAv+Byv-8=0, (212) 
onde 
Vv= ( ê ) ER, 
y 
5 —2 
1=(02 8) 
e 


B=(4V2 1242). 


A matriz A é a matriz simétrica que representa a forma quadrática 
a(x, y) = 2x2 + 4xy +2yº com respeito à base canônica. Não é difícil ver 
(exercício para o aluno) que os autovalores da matriz A são A, =4e A =0,e 
os autovetores normalizados são 


u-(1a 


autovetor associado ao autovalor Ay = 4 
1/v2 ) gs 
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wo —( Riso 


autovetor associado ao autovalor À, = 0. 
1/2 ) , 2, 


Como (u;, u2) forma uma base ortonormal de R?, então 


ia si PINO pa 
p=Iu u)= (15 1/v2 : 


é uma matriz ortogonal que diagonaliza a matriz 4, e a matriz diagonal corres- 


pondente será 
40 
pelc o 


A forma diagonal de q é dada por 


E, também, D = P'AP. 


4 0 x 
g(x1, 71) = (1 71) 0 ) E ) 
= di 
isto é, 
vVAv=41, 
onde 
so) 
V= e Vi=— , 
y y1 
com 


v=Py, ou vi=P'v. 


Como det(P) = 1, observe que v = Pv; é uma rotação. A forma linear se 
transforma em 


Bv =B(Pvi) 


= BPv; 


e mia vhB ) (a) 


= (16 9 (5) 


= 16x,+8y1. 


Substituindo 
vAv=4x e By= 16x +8y 
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em (27.2), obtemos 
4 + 16x) +87, —8=0, 


ou, simplificando, 


+40 +2y-2=0. (27.3) 


Completando o quadrado na variável x1, 
Z+m=(n+2)-a. 
E, substituindo em (27.3), obtemos 
(x1+2)2-44+2y/-2=0, 


ou 
(1 +2/+2(1-3)=0. (27.4) 


Essa equação já é uma forma bem mais simples da cônica inicial e já se pode 
identificar a equação de uma parábola, mas ela ainda pode ser mais simplifi- 
cada. Realizando a mudança de variáveis em (27.4) dada por 


À w =x +2 
y»2=)1—3, 


que representa uma translação no R?, obtemos 
2 
Xy = —2y 2» 


que representa a cônica inicial aos novos eixos-x2y2. Nessa forma, identifi- 
camos facilmente a equação de uma parábola, como ilustra a Figura 27.5. 


XY 


X, 


Figura 27.5: A parábola x = —2y». 
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PROCEDIMENTO PARA SIMPLIFICAR A EQUAÇÃO DE 
UMA CÔNICA 


Seja a cônica I” dada pela equação 
a? +bay+cy +dxrey+f=0. 
Podemos reescrevê-la na forma matricial, 


vVAv+Bv+f=0, 


onde 
g(x,v) =ax +by+o” 
= a b/2 x 
ela O) 
=v'Av, 
e 
f(x,y) =dx+tey 
X 
= (d 
(a (5) 
— By, 
com si2 
a 
Ei e c : 
B =(d e) 
e 


A ideia principal do procedimento a seguir consiste em realizar uma 
rotação nos eixos-xy, de modo a eliminar o termo cruzado bxy. 


1º Passo: Encontrar uma matriz ortogonal P = [u, u2] que dia- 
gonalize A. Lembre que as colunas de P formam uma base (uy, u2) 
ortonormal de autovetores da matriz À para o R2. Assim, 


M O 
— pt = 
D=PAP com D=( 0 o: 


onde À e À» são os autovalores da matriz A associados aos autovetores 


u, € u,, respectivamente. 


2º Passo: Permutar as colunas de P, caso seja necessário, de modo 
que se tenha det(P) = 1. Isso garante que a transformação ortogonal 


x 
v=—Pv;,, com vi= : ; 
1 


seja uma rotação no plano. 


3º Passo: Obter a equação que representa a cônica I” no novo sis- 
tema de eixos-x,y2. Para isso, observe que 


ax +bxy+cy? =vAv 
= (Pvi)'A(Pvi); onde y = Pv; 
=vi(PAP)vi 
= vi/Dvi 


= (x1 mo (5 des 


= Axf + Aoy7 ; 


dxt+ey =Bv 

= B(Pv1); onde v = Pv; 

= (BP)vi; onde BP=(d, e1) 

41 
= (di e 
(dh 1) ( y1 ) 

= dx +ey. 

Assim, a equação v'Av + Bv + f = O se transforma em 
Ax+oyjtdxs Hen +f=o, 


que é uma equação que representa a cônica I” e não contém termos 
cruzados (em xy). 


Vamos fazer uma breve análise dessa equação. 


1. Considere o caso em que os autovalores são não-nulos: 
M £ 0,M» £ 0. Neste caso, podemos completar os quadrados 
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nas variáveis x1 e y1, obtendo 
Axf +Aoy7 +dix+teyi+f= 
= (Mxj+dixi) + (Aoyf + ei) 


=x +) +F, 


com F € R. Assim, a equação 


Adtiy+dx Hey +f=0 
é transformada em 
A + Ao) +F =0. 
Note que 


a. Se AM > 0,1 > 0, então a cônica T” será uma elipse, caso 
F <0; ou um ponto ((x2, y2) = (0, 0)), caso F = 0; ou o 
conjunto vazio, caso F > 0. 


b. Se A <0,Ap <0, então a cônica T' será uma elipse, caso 
F >0; ou um ponto ((x2, y2) = (0, 0)), casoF =0; ou o 
conjunto vazio, caso F < 0. 


c. Se Ay <0< 2,, então a cônica T' será uma hipérbole, caso 
F 0; ou um par de retas concorrentes, caso F = 0. 


2. Considere o caso de um autovalor nulo, digamos, À = 0e À» £0 
(necessariamente À» < 0). Novamente, completando o quadrado 
na variável y1, obtemos 


doyitdmn+eyn+f =(Ayjtey)+dm+f 
= Aoy5 +dyyo+F. 
Assim, a equação inicial da cônica T” fica transformada em 
My +dx+F=0. 
Note que 


a. Sed; 0, então TI será uma parábola. 


b. Se d; = 0, então T” será um par de retas paralelas, caso 
Ao -F <0; ou uma única reta, caso F = 0; ou o conjunto 
vazio, caso À, -F > 0. 
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3. Ocaso À = 0e À £ O é análogo ao anterior. 


E importante observar que nunca poderemos ter À; = À» = 0, pois 
estamos supondo que a forma quadrática associada é não-nula. 


Veja, também, que 


Md GR 0 


0) 


2 
Portanto, À, - À» tem o mesmo sinal de ac — 7º que por sua 


vez tem o mesmo sinal de 4ac — b?. Assim, podemos refazer a 
análise anterior em função do discriminante b? — 4ac da forma 
quadrática. 


Teorema 27.1. 

Dada a cônica de equação ax? +bxy +cy? +dx+ey+ f =0, então 
esta cônica representa: 

a. uma elipse, um ponto ou o conjunto vazio, caso 


b? —-4ac <0; 


b. uma parábola, duas retas paralelas ou uma única reta, caso 


b? -4ac =0; 


c. uma hipérbole ou duas retas concorrentes, caso 


b? — 4ac > 0. 
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Autoavaliação 


Esta aula constitui uma excelente aplicação dos conceitos vistos 
nas aulas anteriores. No entanto, pressupomos que você tenha al- 
guns conhecimentos acerca das equações de cônicas tradicionais, 
como elipses, parábolas e hipérboles. Conhecendo essas equações 
e com o conhecimento adquirido das últimas aulas, você não deve 
encontrar muita dificuldade para compreender os conceitos apre- 


sentados aqui. No entanto, como esta aula reúne muitos conheci- 
mentos matemáticos, você deve ser persistente na leitura dos exem- 
plos e do procedimento apresentado, sempre recorrendo ao tutor no 
caso de encontrar uma dificuldade maior. Na próxima aula, tratare- 
mos de equações semelhantes, agora com três variáveis ao invés de 
duas, mas o procedimento será exatamente o mesmo, ou seja, diago- 
nalizar uma forma quadrática e completar quadrados até simplificar 
a equação ao máximo. 


Exercício 27.1. 


1. Dada a cônica de equação 2x? — 4xy—y? — 4x — 8y+14=0, aplique 
o procedimento apresentado nesta aula, simplificando a equação 
ao máximo e identificando a cônica apresentada. 


Aula 2 S 


QUÁDRICAS 


Ao final desta aula, você deverá ser capaz de: 


compreender o conceito generalizado de uma quádri- 
ca; 

aplicar os conceitos apresentados em casos particu- 
lares. 


Álgebra Linear II | Quádricas 


Pré-requisitos 
Aulas 22, 25,26 e 
27. 
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QUÁDRICAS 


Esta aula é uma continuação da aula anterior sobre cônicas; nela 
estudaremos as superfícies quádricas no espaço Rº. Mais precisamente, 
vamos estudar alguns conjuntos de Rº cujas coordenadas, com respeito 
à base canônica, satisfazem uma equação do tipo 


a? +by+cê+dxytexz+fyztexthy+kz+p=0. 


Usando novamente os resultados sobre diagonalização de formas 
quadráticas, iremos simplificar essa equação e descrever as superfícies 
mais simples que ela pode representar. 


Definição 28.1. 


Uma superfície quádrica, ou, simplesmente, uma quádrica, é o con- 
junto de pontos de Rº cujas coordenadas (x, y, z) satisfazem uma 
equação da forma 


a +by e? +dxytexz+ fyztgxthytiz+p=0, (28.1) 


onde os coeficientes a, b, c,..., k, p são números reais e pelo menos 
um dos coeficientes a, b, c, d, e, f é não-nulo. 


Observe que a equação (28.1) contém uma forma quadrática não- 
nula em Rº, 


ELO PRA Rae a +by +cr +dxy+exz+ fyz, 


uma forma linear em Rº, 


tr, 2)=gx+hy+kz, 
e o termo constante p. 


Apresentaremos a seguir os exemplos mais comuns de superfícies 
quádricas. 


> 


c. Hiperbolóide de duas folhas 


x y 7 E 
CRC A 


b. Hiperbolóide de uma folha 


d. Cone elíptico 


2 2 2 
(=+a-5-0). 
a b2 cd 


y 
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e. Parabolóide elíptico 


(a,0,0) 


g. Cilindro elíptico 


x y 
(= += : 


f. Parabolóide hiperbólico 


h. Cilindro parabólico 


Figura 28.1: Gráficos de Quádricas. 


Observe que a equação (28.1) também pode representar um con- 
junto vazio (por exemplo, x? +y? + 1 = 0), um único ponto (por exem- 
plo, x2+y2 + (z— 1)? = 0), um plano (por exemplo, z? = 0), dois planos 
paralelos (por exemplo, z? = 4) ou dois planos secantes (por exemplo, 
xz = 0). Nestes casos, as quádricas são ditas degeneradas. 


Assim como foi feito para as cônicas, mostraremos que através de 
uma mudança de coordenadas podemos reduzir a equação (28.1) de 
modo que a quádrica seja identificada como sendo de um dos tipos des- 
critos. Esse problema é o de classificar a quádrica. 


Sempre que a quádrica for representada por uma equação que não 
contém termos em xy, xz, YZ, X, ) € z, dizemos que a equação está na 
forma canônica e que a quádrica está na posição canônica. A presença 
de termos cruzados da forma xy, xz ou yz na equação (28.1) indica que 
a quádrica sofreu uma rotação com respeito à posição canônica, e a 
presença de termos da forma x, y ou z indica que a quádrica sofreu uma 
translação com respeito à posição canônica. 


Como foi feito no caso das cônicas, vamos desenvolver um pro- 
cedimento para representar uma quádrica na forma canônica. A ideia 
principal do procedimento consiste em obter um novo sistema de co- 
ordenadas x1y1Z1 de modo que não apareçam os termos cruzados x,y1, 
X1Z1 C)IZI- 


Vamos, primeiramente, expressar a equação (28.1) na forma matri- 
cial. Temos, 


g(x,7,2) =a2+by +crl+dyy+exz+ fyz 


a d/2 e/2 x 
=» 0)| d2 b ff y 
EU DO ip do ci Z 


onde 


s 
Q 
a, 

rd 
9) 
Q 

a 
9) 


Observe também que 
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t(x,7,2) =gx+Hhy+kz 


x 
=(g hk|y 
Z 
= By, 
onde 
Do fish) 


Substituindo g(x, y, 2) = v'Ave (x, y, z) = Bv em (28.1), obtemos a 
forma vetorial da quádrica, 


vVAv+4Bv+p=0. (28.2) 


PROCEDIMENTO PARA SIMPLIFICAR A EQUAÇÃO DE 
UMA QUADRÁTICA 


Seja I” a quádrica representada pela equação (28.1), 
at+by+cê+dxytexz+fyztexthy+kz+p=0, 
cuja forma vetorial é a equação (28.2), 
vVAv+By+p=0. 
1º Passo: Encontrar uma matriz ortogonal P = [u, uz us| que di- 
agonaliza A. Como já foi visto várias vezes ao longo do curso, lembre 


que as colunas de P formam uma base ortonormal (u,, u>, u3) de au- 
tovetores da matriz A para o Rº. Assim, 


M O O 
D=PAP com D= | 0 Ah O |, 
O O ds 


onde À1, À» e À3 são os autovalores da matriz A associados aos 
autovetores u1, U2 e u2, respectivamente. 


2º Passo: Permutar as colunas de P, caso seja necessário, de modo 
que se tenha det(P) = 1. Isso garante que a transformação ortogonal 


v=—Pv;,, com vi= | y1 


seja uma rotação no plano. 


3º Passo: Obter a equação que representa a quádrica TI” no novo 
sistema de eixos x1)1Z1. Para isso, observe que 


ax +by +c7? +dxy+exz+ fyz=V'Av 
= (Pvi)'A(Pvi); onde v = Pv; 


=vi(PAP)vi 
=v/Dvi 
M 0 0 X1 
= 7 2)[ 0 d& 0 y1 
0 0 Às Fál 
E Maxi + Aoy7 + As 


gx+hy-+kz = Bv 
= B(Pv1); onde v = Pvi 
= (BP)vj; ondeBP=(g; hy 1) 


=(gg h k)| 9d 


=gima+Hhyi+ka. 


Assim, a equação 
vAv+Byv+p=0 


se transforma em 
AMZ + Ao + As +gm +hyyi+kiga +p=0. 


Essa equação representa a quádrica I” e não contém os termos cruzados 
X1Y1, X1Z1 CVIZI. 


4º Passo: Completando os quadrados em x1, y1 € z1, obtemos 
(Mt + g1x1) + (Aoy7 + vi) + (Ay + +Hhizi) +p=0 
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M 


h k 
MOS + 21) ++ Ton) +da(ed + quai) + p= 0 
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h4 
M 1614590) + (1 oa 


Passando para as novas variáveis 


us o Mo 
Ap" Z2=2 + ; 


X2 =X + = y1+ 


o 5 2 = 


obtemos a equação 
As + Ao + As25 +pi =0. 


Essa equação representa a quádrica I” e não contém os termos cruzados 
X2)2, X222 € Y2Z2 nem os termos em x», y2 e zo. Portanto, é uma equação 
na forma canônica. 


Exemplo 28.1. 


Descreva a superfície quádrica cuja equação é dada por 


4 +47) +42 +4xy+4x244y2—3=0. 


Solução: 


Reescrevendo essa equação na forma matricial, temos 


vAv-3=0, (28.3) 
onde 
x 422 
v=| y eA=|[2 42 
Z DD A 


Deixamos para você o exercício de calcular os autovalores e os autovetores 
correspondentes da matriz A. Obtemos: 


e À = 2: é um autovalor com multiplicidade algébrica 2 e autovetores 


associados 
=1//2 —1/v6 
u= 1/2 | em=| -1//6 |; 
0 2/v'6 


e À = 8: é um autovalor com multiplicidade algébrica 1 e autovalor as- 
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sociado 


1/V3 
u=| 1/3 
1/V'3 


Como (u,, u2, us! forma uma base ortonormal de Rê, temos que 
3 


1/2 -1/V/6 1/V3 
P=[lu; Us, us] = 1/V2 —1/v6 1//3 
0 2//6 1/3 


é uma matriz ortogonal que diagonaliza a matriz A e a matriz diagonal corres- 
pondente será 


vo) 

I 
SON 
Ono 
o oo 


Vale também que D = P'AP. 


Observe que det(P) = 1, logo P representa uma rotação em Rº. Con- 
siderando 


X X1 
V=|) evi=)| 01 
Xe Z 


e substituindo v — Pv; em v'Av, obtemos 
vAv = (Pvi)'A (Pvi) 
= vi (P'AP)vi 


—v'Dv, onde P'AP — D 


200 X1 
== (x1 »1 z1) 020 »1 

0 08 Z 
=x +23 + 877. 


Portanto, substituindo 
vVAv=24 +27 +84 
na equação (28.3), obtemos 
2% +2yj +87 =3, 


ou, equivalentemente, 
2 2 2 
M x 4 


3/2 3/2 '3/8 
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Observe que essa equação não contém os termos cruzados x1Y1, X1Z1 CY1Z1 
nem os termos em x4, yj € z1. Portanto, é uma equação na forma canônica. 
Identificamos, facilmente, que essa equação representa um elipsóide, como 
ilustra a Figura 28.1.a. 


Exemplo 28.2. 


Identifique a superfície quádrica cuja equação é dada por 


x +2yz— 29422-101 =0. 


Solução: 


Inicialmente, observe que a presença do termo cruzado yz nos levará a 
realizar uma rotação de eixos, e a presença dos termos lineares z e y, a realizar 
uma translação de eixos. 


Reescrevendo essa equação na forma matricial, temos 


vAv+Bv-101=0, (28.4) 
onde 
x —1 0 0 
v=|y |, 4=[ 001 |[eB=(0 -v2 42). 
Z 010 


Deixamos para você, novamente, o exercício de calcular os autovalores e 
os autovetores correspondentes da matriz A. Obtemos: 


e À = —1: autovalor com multiplicidade algébrica 2 e autovetores asso- 
ciados 
1 0 
u=|0 e uw= 1/ v2 : 
0 —1/v2 


e À, = 1: autovalor com multiplicidade algébrica 1 e autovalor associado 


0 
us =| 1/2 
1//2 


Como (u,, u2, us) forma uma base ortonormal de Rº, então 
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1 0 0 
P=[u; us, us] = 0 1/42 1//2 
0 -1//2 1/42 


é uma matriz ortogonal que diagonaliza a matriz A e a matriz diagonal corres- 
pondente será 


-1 00 
D= 0 —1 0 
0 01 
Vale também que D = P'AP. 


Como no Exemplo 28.1, det(P) = 1, logo P representa uma rotação em 
Rº. Considerando 


X X1 
E O evi=| ) 
Z Z 


e substituindo v — Pv; em v'Av, obtemos 
vVAv = (Pvi)'A (Pvi) 
= vi (P'AP)vi 


— v Dvi, onde P'AP — D 


—1 0 0 X1 
=(x 121) O =) y1 
0 1 2 
=-A ota, 
e, substituindo v = Pv, em Bv, obtemos 
Bv = B(Pvi) 
= BPv; 
1 0 0 X1 
= (0 -v2 vD)| 0 1/N2 1/2 Yi 
0 —1/v2 1/v2 z1 
x 
= (0-2 0)| 9 
z1 
— —21. 


Portanto, substituindo 


vAv=-"-y+ã e Bv=-2y 


em (28.4), obtemos 
= ye 2yr = 101 
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Agora, completando o quadrado na variável y,, temos 
-4 + —(4+2n)=101, 


o que nos dá 


n—[01+12-1] +2j= 101, 


e, portanto, 


ou, equivalentemente, 


2 2 2 
M (1 +1) Ei 

= = ECO Soo. 28. 
102 102 o 102 Cas 


Essa equação já é uma forma canônica para a quádrica inicial e já se pode 
identificar a equação de um hiperbolóide de duas folhas, mas ela ainda pode 
ser mais simplificada. Realizando a mudança de variáveis dada por 


X2 — X1 
32 = 151 
221 


que representa uma translação no Rº, a equação (28.5) se transforma em 


2 
Z2 
Lt =] 
+10 j 


ERA 
102 102 


que representa a quádrica inicial aos novos eixos x2y2z2. Nessa forma, identi- 
ficamos novamente a equação de um hiperbolóide de duas folhas, como ilustra 
a Figura 28.1.c. 


Autoavaliação 


Terminamos o estudo das cônicas em R? e das quádricas em Rº, 
que constituem uma excelente aplicação da diagonalização das for- 


mas quadráticas. É importante que você reveja o procedimento de 
simplificação dessas equações e compreenda os cálculos realizados 
nos exemplos. Também é importante que fique clara a interpretação 
geométrica de cada mudança de variáveis realizada. 
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Exercício 28.1. 


Obtenha uma forma canônica de cada quádrica abaixo e identifique 
a quádrica. 

1. 2xy—-4/2x+42/2y+7-9-0. 

2. 2xy+2x2+4+2y2—6x—6y—-472-9=0. 


3. 1x2 + 772 4 107º — 2xy — 4x7 4 4yz — 12x + 1274 607+ 66 = 0. 
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Aula 2 Õ 


AUTOVALORES COMPLEXOS 


Ao final desta aula, você deverá ser capaz de: 


compreender o conceito de autovalor complexo; 


aplicar os conceitos apresentados em casos particu- 
lares. 


Álgebra Linear II | Autovalores Complexos 


Pré-requisitos 
Aulas 3 e 5. 
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AUTOVALORES COMPLEXOS 


Vimos logo na Aula 3 que, dada uma matriz A E M, (IR), seu polinô- 
mio característico p(x) é um polinômio de grau n com coeficientes reais 
e, portanto, possui um total de n raízes, contando suas multiplicidades 
e as raízes complexas. Nesta aula, estudaremos alguns exemplos de 
matrizes reais com autovalores complexos. 


Inicialmente, vamos relembrar alguns conceitos sobre números com- 
plexos. Denotamos o conjunto dos números complexos por C e repre- 
sentamos por 


C=(a+bila,beRe i=vTh 


A igualdade de números complexos é definida por 
a+bi=c+diseesomentesea=ceb=d. 


A adição e a multiplicação de números complexos são definidas por: 


a. (a+bi)+(c+di) = (a+c)+(b+d)i; 


b. (a+bi)-(c+di) = (ac — bd) + (ad + bc) i, 


para todos os a, b, c, dE R. E fácil verificar que todas as propriedades 
de corpo dos números reais continuam válidas para os números com- 
plexos. 


Definimos o conjugado de um número complexo z = a + bi como 
sendo o número complexo 7 = a — bi. 


A teoria de espaços vetoriais e de álgebra matricial desenvolvida 
no caso de componentes reais e escalares reais se aplica também para 
componentes e escalares complexos. Por exemplo, o espaço vetorial C? 
é definido por 

C=f(z,w|zw el 


com as operações usuais 


a. (21, wi) + (22, w2) = (21 +22, W1 + W9); 


bravo (tra! 


onde Z, Z1, W1, 22; W2 E (os 


Assim, dada uma matriz A E M,(C), um número complexo À e € 
é um autovalor (complexo) da matriz A se existe um vetor não-nulo 
veC” tal que 
Av= Ay. 


Dizemos que v é um autovetor (complexo) associado ao autovalor 


AEC. 


Exemplo 29.1. 


Discuta a diagonalização da matriz 
0 —1 
Ai ( a ) | 


Solução: 


Sabemos, do nosso estudo de rotações no plano, que essa matriz correspon- 
de a uma rotação de 7/2 radianos no sentido anti-horário em torno da origem 
do plano cartesiano R?. Assim, fica claro que nenhum vetor não-nulo v € R? é 
transformado, pela ação da matriz A, num múltiplo dele mesmo. Assim, a ma- 
triz A não possui autovetores em RR? e, consequentemente, não tem autovalores 
reais. De fato, o polinômio característico de A é 


p(x) =det(xb — A) 
x 1 
—l x 
El: 


Esse polinômio só possui as raízes complexas À = ie o = —i. 


No entanto, considerando A com matriz complexa, isto é, AE M>(C), 
hM =ie A = —i são autovalores complexos da matriz A, pois os vetores 
v=(1,-i),v=(1,i) € C2, e satisfazem 


me (PDM) 
m((DA()-m 


Assim, vi = (1,—i) é um autovetor associado ao autovalor Ay = i, e 
v2 = (1, i) é um autovetor associado ao autovalor À, = —i. 


Como a matriz não possui autovalores reais, ela não é diagonalizável en- 
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quanto matriz real. No entanto, como ela possui dois autovalores complexos 
distintos, a matriz A é diagonalizável quando considerada como matriz com- 
plexa. Mais ainda, considerando as matrizes P, D e M>(C) dadas por 


temos 
PDP! = y (o Boo 0] 


(10) (1 dê) 


Pd 


isto é, A = PDP”!. Portanto, no caso complexo, a matriz A é semelhante à 
matriz diagonal D. 


Exemplo 29.2. 


Dada a matriz 
E 0,5 —0,6 
TE is UE a 


determine os autovalores de A e uma base para cada autoespaço. 


Solução: 


Obtendo o polinômio característico da matriz A, 


p(x) =det(xb — A) 
x—0,5 0,6 
—0,75 x—1,1 
=(x—0,5)(x— 1,1) — (0,6)(—0,75) 


=x —-1,6x+1. 


Calculando as raízes desse polinômio quadrático, obtemos 


M =0,8-0,6ie A =0,8+0,6i. 
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Considerando o autovalor À, = 0,8 -— 0, 61, queremos obter y = (z, w) E C? 
não-nulo tal que 
Av= My, 


0,5 —0,6 Z = = R Z 
( us did ) ( ay HU 0) ( w ). 


o que nos dá o sistema linear 


ou seja, 


(-0,3+0,6i)2-0,6w =0 
0,752+(0,3+0,6i)w = 0. 


Como os autovalores são distintos, cada autoespaço tem dimensão 1; portanto, 


as equações do sistema anterior são dependentes. Assim, basta considerar uma 
das equações; por exemplo, da segunda equação, temos 


z=(-0,4-0,8i)w. 


Escolhendo w = 5 (para eliminar a parte decimal), obtemos z = —2 — 4i. As- 
sim, uma base para o autoespaço associado ao autovalor Ay = 0,8 — 0,61 é 


dada pelo vetor 
E —2-—4i 
| 5 cá 


Analogamente, para o autovalor À, = 0,8 +0,6i, obtemos o autovetor 


fura 
Mor= 5 , 
As es DO 06 =D di 
gos Jd 5 
= é dA 
FÃ dE 


= (0,8+0,6i) ( 


pois 


-244i 
5 


E Aovo R 


Observe que a matriz A é semelhante à matriz diagonal 


p=(% 0) (08-06 0 
ar gs ar 0 0,8+0,6i )" 
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Autoavaliação 


Não é nosso objetivo generalizar toda a teoria de diagonalização de 
matrizes reais para o caso complexo; apesar disso, desejamos pro- 
porcionar novas e importantes aplicações da Algebra Linear. Muitos 


problemas envolvendo matrizes com autovalores complexos apare- 
cem naturalmente em Engenharia Elétrica, em Física e na área de 
Sistemas Dinâmicos de um modo geral. Essa discussão costuma ser 
feita num curso avançado de Álgebra Linear. Portanto, nosso obje- 
tivo foi apenas o de apresentar a você alguns exemplos elementares. 


Exercício 29.1. 


1. Determine os autovalores e uma base para cada auto-espaço da 


matriz 
1 —2 
(e) 
2. Calcule os autovalores e autovetores da matriz 
a —b 
do a) 
ondea, beR coma £0o0ub 0. 


3. Dada a matriz A E M,(R) com autovalor À E C, mostre que À 
também é autovalor da matriz A. 
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EXERCÍCIOS RESOLVIDOS — 3º PARTE 


Nas próximas aulas apresentaremos uma série de exercícios resolvi- 
dos sobre a segunda parte do curso. Esses exercícios o ajudarão a con- 
solidar os conceitos apresentados nas aulas anteriores. 


A nossa orientação é que você primeiro tente resolver cada um dos 
exercícios, usando, se necessário, as anotações das aulas anteriores, e, 
só depois de obtida a sua própria solução, compará-la com a solução 
apresentada aqui. Caso você não consiga resolver algum exercício, não 
se aflija, leia atentamente a solução correspondente. Se você ainda tiver 
dificuldade, não hesite em procurar ajuda de seu tutor. 


Exercício 30.1. 


1. Determine a matriz, com respeito à base canônica, da projeção 
ortogonal sobre a reta y = x. 


2. Determine as projeções ortogonais dos pontos P; = (1,0, 1) e 
P =(1,1,1)sobreo planox+y-—z=0. 


3. Determine o valor das constantes a, b, c,d E R para que 


l a+b b 3 b=€ 2d +3 
A=|2 0 4 cB=| 3 5 l 
3 4 3 d b+c 0 


sejam matrizes simétricas. 


4. Dadas as matrizes simétricas 4, B € M,(R), mostre que AB + BA 
também é uma matriz simétrica. 


5. Dadas as matrizes 4, BE M, (IR) tal que A é uma matriz simétrica, 
verifique que B'A B é uma matriz simétrica. 


6. Dados a, b E R, com b £ 0, encontre uma matriz ortogonal P que 


diagonaliza a matriz A = E ; ) , isto é, tal que D = P'AP. 

7. Seja T: Rê — Rº um operador autoadjunto com autovalores asso- 
ciados À = 3 e À» = 4; suponha que v, = (1, 1,1) e 
v2= (2,0, 1) são dois autovetores associados ao autovalor 
A = 3. Determine um autovetor associado ao autovalor À, = 4 e 
uma base ortonormal de autovetores de T. 


8. Para cada matriz abaixo, determine uma matriz ortogonal P e uma 
matriz diagonal D tais que A = PDP”. 


Ed cid 31000 
E SO 13000 
a A- b4=/1 00211 
0000 
0-0 00121 
001172 
SOLUCÕES 


1. Denotamos por T:R? — R? a projeção ortogonal sobre a reta 
y =x, como ilustra a Figura 30.1. 


Figura 30.1: A projeção ortogonal sobre a reta y = x e a base ortonormal B. 


Vamos primeiro determinar uma matriz que representa T com res- 
peito a uma base ortonormal 5 = (u,, u>+. Sejam: 


u, = (1/V2, 1/2) vetor unitário paralelo à reta y = x; € 
w =(—1/v2, 1/2) um vetor unitário normal à reta y = x. 


Como 
T(ui) =u,=1l-u,4+0-u 


T(u,) =0=-0-u,+0-u, 
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temos que 


s=Ino=(0 a 


Assim, a matriz A que representa 7 com respeito à base canônica 
é dada por 
A=PBP, 


onde 
- = 1/2 —1/v2 
fes ul (a o 


Como P é uma matriz ortogonal, temos que 


pi=p=( Ez a): 


= /y2 1/2 
portanto, 
A=PBP| 
(ia ia) (a 6) (a a) 


(ih 1a) 


2. Seja T:Rº -—» Rº a projeção ortogonal sobre o plano 


m:x+y—z= 0; precisamos determinar a matriz A que repre- 
senta essa projeção com respeito à base canônica. Novamente, 
vamos primeiro obter a matriz que representa T com respeito a 
uma base ortonormal 8 = (u,, uz, us). Veja a Figura 30.2 


Figura 30.2: Uma base ortonormal B. 


Considere os seguintes vetores: 
u, = (1/2,0, 1/2) um vetor unitário paralelo ao plano 7, 
w=(-1//6,2/v6, 1/6) um vetor unitário ortogonal au e 
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paralelo ao plano 7 e 
us = (1/43, 1/43,-1/3) um vetor unitário normal ao plano 


7. 
Como 
T(u) =uy =1-u,+0-w +0-us; 
T(u)=u, =0-u, +1-u+0-us; 
e 
T(us)=0=0-u, +0-u,+0-us, 
temos que 


100 
B=[T]p=[ 010 
000 


Assim, a matriz A que representa T com respeito à base canônica 
é dada por 


A=PBP, 
onde 
1//2 —1/V6 1/43 
P=[u; us, mi 0 2/v'6 ral 
1/2 1/6 -1/V3 


Como P é uma matriz ortogonal, temos que 


1/2 0 1/v2 
remo (ING 2/6 nã) 
1/43 1/43 —1/43 
portanto 
A=PBP 
1/2 -1//6 1/V3 100 
0 2/6 na) (010). 
1/2 1/6 —1/V3 000 
1/2 0 1/v2 
(rá 2/V'6 na 
1//3 1/3 —1//3 
2/3 —1/3 1/3 
=| -1/3 2/3 1/3 |. 
1/3 1/3 (5) 
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Assim, as imagens dos pontos P, e P,, sob a ação da projeção 
ortogonal sobre o plano 7x, são obtidas por multiplicação de ma- 


trizes: 
2/3 —1/3 1/3 1 1 
AP, = —1/3 2/3 1/3 O 1=|1 0 |: 
1/3 1/3 23 1 1 
2/30 AS dB 1 2/3 
AB = —1/3 2h3 M/S 1 E DAS 
1/3 1/3 2/3 1 4/3 


Portanto, temos que a projeção de P, é P, e a projeção de P, é 
(2/3, 2/3, 4/3). 


3. Lembre que uma matriz A é simétrica se e somente se 4 = A'. 
Assim, para a matriz 


1 atb b 
AS) D SO A a 
3 A 3 
temos A = A! se e somente se 
l a+b b 
o, 0 4 = a+b 0 4 |, 
o ct b 43 
ou seja, se e somente sea+b =2 eb =-3,ou, ainda, a = —1 e 
b=3. 
Para a matriz 
5 b—-c 2d43 
B=[3 5 1 
d b+c (0) 
temos B = B' se e somente se 
5 b—-c 2d+43 3 3 d 
o O 1 = b-c 5 b4c |, 
d b+c 0) 2d4+3 1 0 


ou seja, see somenteseb—-c=3, b+c=le 2d+3=d,ou, 
ainda,b=2,c=-le d=-3. 
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4. Sendo A e B matrizes simétricas, temos A = A' e B= B”. Portanto, 


(AB+BA)' = (AB) +(BA) 
— B'A! + A'B' 
= BA +AB 
= AB+4-BA. 


Portanto, a AB + BA também é uma matriz simétrica. 
. De fato, temos que 


(B'AB)' e B'A'(B')' 
— B'AB; 


logo, B'AB também é uma matriz simétrica. 


. Como A é matriz simétrica, existe uma matriz ortogonal P que 
diagonaliza a matriz A. Lembre que as colunas de P são autove- 
tores unitários da matriz 4. Portanto, precisamos calcular os au- 
tovalores e os respectivos autovetores da matriz A. Seu polinômio 
característico é dado por 


p(x) =det(xhb — A) 


| d=4 —b 
—b x-a 
= (ra)? -(-b) 


=x -2ax+(a2 —b?). 


Portanto, os autovalores são A = atbe Ap =a-—b. Como 
b 0, segue que A * À». Vamos, agora, ao cálculo dos autove- 
tores. O autovetor associado ao autovalor À = ab é um vetor 
uy = (x,y) E Rº que satisfaz 


(Md —A) uy, — 0, 


b bNV/2) (0 

E RO DADA 
Como b + 0, obtemos x = y. Assim, uma escolha de uy = (x, y) 
que seja vetor unitário é dada por uy = (1/v2, 1/v2). Como 


M *ho e a matriz A é simétrica, então todo autovetor 
u, = (x, y) € R? associado ao autovalor À = a — b é ortogonal 


ou seja, 
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ao vetor uj. Portanto, podemos escolher us = (—1/v2, 1/v2). 
Assim, a matriz 


«VAND SInO 
P=Im vl= (15 108) 


diagonaliza a matriz A, isto é, 
D=P'AP 


(oa na) (6 a) Cipa pa) 
=(05" e 


é uma matriz diagonal semelhante à matriz A. 


. Seja V3 E R? um autovetor associado ao autovalor À = 4. Como 


T é um operador autoadjunto e os vetores vj e v2 são linearmente 
independentes, devemos ter v3 ortogonal a vj e v>. Como esta- 
mos em Rê, v3 é paralelo ao vetor v; x v2: portanto, podemos 
considerar 


Vv3=Vixv=(1,1,-2). 


Observe que para os autovetores vj e v> associados ao autovalor 
M = 3 temos 


(Vi, V2) = 1EM 0, lj, (1, de D) 
=D aco op doi 
320; 


logo, vj e v> não são ortogonais entre si. Para construir uma base 
ortogonal de autovetores, consideramos os vetores v2, v3 e um 
novo vetor w, com w ortogonal a v> e v3, por exemplo, 


w=vxv=(—1,5,2). 


Normalizando esses vetores, obtemos uma base ortonormal de 


autovetores 8 = (u,, uz, us), dada por: 
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Ra (5 1 =) 
2 Tysl 6 V6' V6)' 


a. Sendo 


O ma O 
Ow mm 


0000 


seu polinômio característico é dado por 


p(x) =det(xly — A) 


$=B =1 00 
o EM 300006 
E 

0 0 0x 


=x (2 — 6x8) 
=x (x—D)(x=4). 
Logo, seus autovalores são: 
e À — 0, com multiplicidade algébrica 2; 


e À» =2, com multiplicidade algébrica 1; e 
e À3 =4, com multiplicidade algébrica 1. 
Vamos, agora, calcular uma base ortonormal de autovetores 


de A. Para o autovalor À = 0, sabemos que os autovetores 
associados v = (x, y, Z, t) E Rº satisfazem 


(0-4 — A)v= 0 
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isto é, satisfazem o sistema linear homogêneo 


8: «sf 040 x 0 
sf «85000 a o E, 
O 000 z|"|o 
O Brad t 0 


Escalonando a matriz associada desse sistema linear, no caso, 
a matriz —A, obtemos as soluções 


x=0, y=0 e z,t arbitrários. 


Portanto, escolhendo ora z = 1,t=0, eoraz=0,t=1, 
obtemos que 


u=(0,0,1,0) w=(0,0,0,1) 


formam uma base ortonormal do autoespaço associado ao 
autovalor À = 0. 


Para o autovalor À, = 2, sabemos que os autovetores associ- 
ados v = (x, y, 2, t) E Rº satisfazem 


(2 7) — A)v — 0, 


isto é, satisfazem o sistema linear homogêneo 


Es, = sb x 0 
ei e fiipdo vi lo 
6: 0200 z| 10 
= ig 0 2 t 0 


Escalonando a matriz associada desse sistema linear, obte- 
mos as soluções 


y=-x e z=t=s=0, com x arbitrário. 


Portanto, escolhendo x = 1/ v2 e, consequentemente, 
j= —1/v2, obtemos que 


1 — 
U3 =| TR 0, 0) 
(a v2 
forma uma base ortonormal do autoespaço associado ao au- 
tovalor À» = 2. Finalmente, para o autovalor À3 = 4, os au- 


tovetores associados v = (x, y, Z, t) E Rº satisfazem 


(4-1 — A)v =— 0, 


ou seja, satisfazem o sistema linear homogêneo 


» 200 


l 0 
—1 0 
0 “10 
0 0 


Escalonando a matriz associada desse sistema linear, obte- 
mos as soluções 


y=x e z=t=0, comx arbitrário. 


Portanto, escolhendo x = 1/ v2 e, consequentemente, 
ve dá 2, obtemos que 


l l 
u= (22:00) 

forma uma base ortonormal do autoespaço associado ao au- 
tovalor Às = 4. Como a matriz A é simétrica, observe que 
os autovetores associados a autovalores distintos são ortogo- 
nais. Assim, B = (uy, u2, us, u4! é uma base ortonormal 
de Rº formada por autovetores da matriz A. Portanto, a ma- 
triz ortogonal P, 


0 0 1/2 1/2 
= “100 1/42 1/42 
P=[u; us us u4] = 10 0 0 , 
01 0 0 
e a matriz diagonal D, 
0000 
0000 
BE | g0 20]? 
0004 


satisfazem A — PDP". 
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b. No caso 
31000 
13000 
A=|0 0211, 
00121 
00112 


seu polinômio característico é dado por 


p(x) =det(xl; — A) 


x—3 —1 O) O) O) 
—1 x-3 0 0 0 
= 0 O x-2 -1 — 
0 O) —) x-2 +— 
O) O) —1 —1 x-2 
ESA qui x—-2 —1 —1 
=p 3 tie 
1 1 x—2 


— (xº —-6x+8)( — 6x2 +9x— 4) 
= (x— Dx—2)(x— 42. 


Logo, os autovalores da matriz A são: 
e À, = 1, com multiplicidade algébrica 2; 
e À») =2, com multiplicidade algébrica 1; e 
e Às = 4, com multiplicidade algébrica 2. 
Vamos, agora, calcular uma base ortonormal de autovetores 


de A. Para o autovalor À, = 1, sabemos que os autovetores 
associados v = (x, y, 2, 1, s) ER? satisfazem 


(1 15 — A)v =— 0, 
isto é, satisfazem o sistema linear homogêneo 


=p es GO 0 0 
== O 0 O 
Do o psi =] 
GQ = sl =] 
E a O | 


Ds nt us 
| 
So 000 


Escalonando a matriz associada desse sistema linear, obte- 
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mos as soluções 
x=0, y=0, z=-t—s com t e s arbitrários. 
Portanto, escolhendo t = 0 e s = —1, obtemos o autovetor 
vi=(0,0,1,0,-1). 


Para obter um segundo autovetor v; = (a, b, c, d, e) asso- 
ciado ao autovalor À; = 1 e que seja ortogonal a v;, devemos 
ter 


a=b=0 
ctd+e=0 
c—-e=0, 


sendo que a última equação segue da condição (vi, v2) = 0. 
Uma solução desse sistema linear é dada por 
v2=(0,0,1,-2, 1). Assim, (vi, v2) é uma base orto- 
gonal do autoespaço associado a À = 1. 
Para o autovalor À, = 2, sabemos que os autovetores associ- 
ados v= (x,y, Z,t,5) € Rº satisfazem 


(2 15 — A)v = 0, 
isto é, satisfazem o sistema linear homogêneo 
-1-1 0 0 0 x 0 
-1-1 0 0 0 y 0 
O 0 0-1 -1 a RE 
O 0-1 0 -1 t 0 
O 0-1 1/0 s 0 


Escalonando a matriz associada desse sistema linear, obte- 
mos as soluções 


y=-—-x e z=t=s=0, com x arbitrário. 
Portanto, escolhendo x = 1, obtemos o autovetor 
v3=(1,-1,0,0,0), 


que forma uma base do autoespaço associado ao autovalor 
hM =2. 

Finalmente, para o autovalor ÀA3 = 4, os autovetores associ- 
ados v= (x, y,2,t,s) ER” satisfazem 


(4-1s — A)v =— 0, 
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ou seja, satisfazem o sistema linear homogêneo 


pa 
o 
o 
o 
Ds nau 
| 
So 000 


Escalonando a matriz associada desse sistema linear, obte- 
mos as soluções 


y=x, s=z e t=z, com x e z arbitrários. 


Agindo como no caso do autovalor À; = 1, obtemos os se- 
guintes autovetores associados ao autovalor À = 4: 


v,=(1,1,0,0,0)evs=(0,0,1,1,1),e eles formam 
uma base ortogonal para o autoespaço associado ao auto- 


valor Ày = 4. 


Assim, (Vi, V2, V3, V4, Vs) é uma base ortogonal de Rº 
formada por autovetores da matriz A. Normalizando os ve- 
tores dessa base, obtemos 


UU, — 


Us — 


U, — 


Us — 


( 
( 
u = (5, 5,0,0,0 
( 
( 


Observe, agora, que 8 = (u;, u5>, u3, u4, us! é uma base 
ortonormal de Rº formada por autovetores da matriz A. Por- 
tanto, a matriz ortogonal P, 


P=[u uw us w us|= 
0 0 1/2 1/2 0 
0 O 1/2 1/02 0 
= 1/v2 1/06 0 O Afy3) 
O -2//6 O 0 1/43 
1/2 1//6 O 0 1/43 


e a matriz diagonal D, 


y 

| 
(os ao HR o 
Sono 
Ee No 
o rooo 
Ss 


0 


satisfazem A = PDP”. Lembre que a ordem dos elementos 
da diagonal principal da matriz D depende da ordem das 
colunas da matriz ortogonal P e vice-versa. 
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Aula 3 1 


EXERCÍCIOS RESOLVIDOS — 4º PARTE 


Ao final desta aula, você deverá ser capaz de: 


aplicar os conceitos e as propriedades vistas nas Aulas 
1729, 


Álgebra Linear II | Exercícios Resolvidos — 4º Parte 


Pré-requisitos 
Aulas 17 a 30. 
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EXERCÍCIOS RESOLVIDOS — 4º PARTE 


Nesta aula, vamos dar continuidade à apresentação de exercícios 
resolvidos sobre a segunda parte do curso. Estes exercícios o ajudarão 
a consolidar os conceitos apresentados nas aulas anteriores. 


Mais uma vez, ressaltamos que você deve primeiro tentar resolver 
cada um dos exercícios, usando, se necessário, as anotações das aulas 
anteriores, e, só depois de obtida a sua própria solução, compará-la com 
a solução apresentada aqui. Caso você não consiga resolver algum e- 
xercício, não se aflija, leia atentamente a solução correspondente e, se 
ainda tiver dificuldade, não hesite em procurar ajuda de seu tutor. Uma 
discussão entre alunos e tutor sobre as soluções encontradas é sempre 
muito proveitosa. 


Exercício 31.1. 


1. Para cada caso abaixo, determine a matriz que representa a forma 
bilinear com respeito à base ordenada especificada. 


a. F:RiÍxRº-R dadaporF(u, v)= (u, v) com respeito à 
base B=(u,um,usk, u=(-2,0, 1),m=(1,2,1)e 
Us — (0, 1,—2). 

b F:RêxR? SR dadaporF(u, v)=(u,a)-(v, b), com 
a, be R”, com respeito à base canônica. 


2. Expresse as formas quadráticas abaixo na forma v'Av, onde a 
matriz A é uma matriz simétrica. 


a. g(x1, x2) = 3x4 4725 


b. g(x1, X2; X3, X4)) = ds + jo +323 +24 +2x/x2 +4x1x3 
+ 6x2x3 +77 x/x4 — 2x9X4 


é gue; do; X3)= El + o — x3 +2x1X2 — 3x1X3 + X2X3 


d. g(x1, x2) = —Tx1x2 
à 
e. g(x1, X2,..,Xn) = (cixy + com + --- + coxn)?, com 
c1, €02,..,Cn ER. 


3. Diagonalize as seguintes formas quadráticas: 


a. q(x, y) = 2xy 
b. g(x,y,2)=2xy+2x2+2yz 


Em cada caso, determine uma matriz ortogonal que diagonaliza a 
forma quadrática. 


. Identifique as cônicas representadas pelas equações abaixo. Em 
cada caso, determine uma matriz ortogonal que diagonaliza a 
forma quadrática. 


a. 2x2 +5y? = 20 

b. x2— 1672 +8x+ 128y = 256 

c. 4x2 —20xy+25yº — 15x—6y=0 

. Identifique as quádricas representadas pelas equações abaixo. Em 


cada caso, determine uma matriz ortogonal que diagonaliza a 
forma quadrática. 


a. 2xy+2xz +2yz— 6x — 6y — 47 = —9 
b. 2xy—6V2x4+10/2y+72-31=0 


6. Seja F a forma bilinear de R? definida por 


F((x1, x2), (y1, )2)) = 2x11 — 3x172 + x2y2. 


a. Determine a matriz A que representa F com respeito à base 
oq =t((1,0), (1,1). 

b. Determine a matriz B que representa F com respeito à base 
B = tias 1, 6! ; RE 

c. Determine a matriz mudança de base P, da base B para a 
base q, e verifique que B = P'AP. 
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1. a. Lembre-se da Aula 25, na qual a matriz que representa a forma 
bilinear com respeito à base B = (uy, u2, us + é dada pela matriz 
A = (a;;), onde a;; = F(u;, u;). Neste caso, temos: 


an = F(u,u,)=(u,,u,))= 
E (2, Ú, D, (2 0, D) = 


az = F(u,u2)= (u,u)= 
= ((—2, 0, D, (1, 2, li =: 


as — F(u, Us) =— (ui, us) = 

= 2 0, D, (0, l; 2 =-2; 
1 = F(u,ui,)=(u,ui)= 

= ((1,2,1),(-2,0,1)=-—1; 
a = F(u,u>)=(u,,u)= 

e Is, 1,(1,2,1)=6; 
da — F(u, Us) = (us, Us) = 

= ((1,2,1),(0,1,-2)=0; 
as = F(us,u,) = (us, ui) = 

= ((0,1,-2),(-2,0,1) = —2; 
az = F(us, u2) = (us, u2) = 

= 40,1,=2),;(1,2,1)=0; 
az = F(us, us) = (us, us) = 


Assim, a matriz A é dada por 


dl, =2 
A=| —1l 6 0 
=P 


Observe que 4 é uma matriz simétrica. 


b. Sejam a = (41, a2) e b = (by, b>) vetores com respeito à base 
canônica. Seja A — (a;;) a matriz que representa a forma bilinear 
F(u, v)=(u,a)-(v, b) com respeito à base canônica. Assim, 
temos: 
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an = F(ere;)=(e,,a)-(ei,b) = 

= ((1,0), (ar, a2)) -((1, 0), (bi, b2)) = abr; 
a» = F(er,e»)=(e,,a)-(e2,b) = 

= ((1,0), (ar, a2)) -((0, 1), (bi, b2)) = aba; 
a = F(e,e)=(e2,a)-(es,b) = 

= (0,1), (ai, a2))-((1, 0), (bi, bo)) = aobi; 


ao =— F(es, e2) = (e2, a) E (eo, b) = 
((0, 1); (ar, a2)) ((0, D, (Dj b>)) = q2ba. 


AULA ES MÓDULO 2 


Portanto, 


ab; ab» 
A= . 
ab; ab» 
Observe que, em geral, a matriz A não é uma matriz simétrica. 


2. Como foi visto na Aula 26, temos: 


a. 
g(x1, x2) = (x1 x2) ( 0 7 ) ( o ) = 3x1 4725 
b. 
g(x1, X2; X3; X4) = 
1 12 o X1 
= (x % ) 1 23 —1 X2 
RAM SD de odio O x3 
TLD, —1 0 1 X4 
=— di +25 = Sã +24 + 2x/x2 + 4x1x3 + 6x02x3 + Txyx4 — 2x2x4 
Cc. 
g(x1, x2, x3) = 
1 1 —3/2 X1 
— (x1 X2 x3) 1 1 1/2 X2 
—3/2 1/2 —1 X3 
=x +25 — x +2x1x2 — 3x1x3 + X2X3 
d. 


q(r, x2) = (x1 22) ( E Es ) ( Es ) = —Tx1x2 
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g(x1, %24.:.5Xn) = (cr cr +=: + cata)” = 


= cixt + cias ++ cixi + 2cicoxixo + 2cic3xx3 RR = 
+2Cn-1CnAn-1Xn 


ct cio ce + ciCn x1 
cio 5 c03 «cn x2 
2 
= (x1 X2 e Xri) C103 0203 03 ver C3Cn X3 
C1Cn Cn C3Ch te: cs Xn 


3. a. Observando os coeficientes de q, vemos que a matriz A que repre- 


senta q na base canônica é dada por 


istio: 
10 

Diagonalizar a forma quadrática q é equivalente a diagonalizar a 
matriz simétrica A. Usando os procedimentos já conhecidos sobre 
diagonalização de matrizes simétricas, os autovalores da matriz A 
são Ay = le Az =-—1. A matriz P será obtida a partir de uma base 
ortonormal de autovetores de 4. Efetuando os cálculos, o que é 
um exercício para você, obtemos 


uia (0 


autovetor associado ao autovalor À = 1,e 
1/N2 ) 


= (UA 


autovetor associado ao autovalor À, = —1. 
1/v2 y 


Como (uy, u>! forma uma base ortonormal de R?, então 


E “(1/2 —/N2 
e ul ( 1a 1/v2 ). 


que representa uma rotação de 7/4 radianos, e a matriz diagonal 
correspondente será 
1 0 
D= ; 
0 —1 


onde D = P'AP. Observe que a forma diagonal de q é dada por 


teo= tm lo 1) (4) 


A 


b. Observando os coeficientes de q, vemos que a matriz A que repre- 


senta q na base canônica é dada por 
011 
A=|[ 101 
110 


Procedendo à diagonalização da matriz simétrica A, deixamos os 
detalhes dos cálculos como um exercício para você, obtemos os 
autovalores A; = —1, com multiplicidade algébrica 2, e À» = 2. 
A matriz mudança de variável P será obtida a partir de uma base 
ortonormal de autovetores de A. Efetuando os cálculos, obtemos 


1/6 
u =| —2/ V6 | autovetor associado ao autovalor À, = —1; 


1/6 
1/2 


= 0 autovetor associado ao autovalor À = —1: 


=1/5/2 
1/3 


us= | 1/3 | autovetor associado ao autovalor À, = 2. 


1/48 


Como (uy, u2, u3 ) forma uma base ortonormal de R?, então 


1//6 1/2 1/43 
P=[u; us, us) = -2/V'6 0 1/3 
1//6 —1/V2 1/43 


é uma matriz ortogonal e a matriz diagonal correspondente será 


-1. 0 0 
D= 0 —-1 0 |, 
0 02 
onde D = P'AP. 
A forma diagonal de q é dada por 
-1. 00 
go), a) =(x Ya 21) 0 —1 0 y1 
0 02 
=--y+d. 


Como P é uma matriz ortogonal e det(P) = 1, então P é uma 
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rotação em Rº. 


. Como a forma quadrática g(x, y) = 2x? + 5y? não contém ter- 


mos em xy, a equação da cônica já está diagonalizada. Podemos 
escrevê-la na forma 


e, daí, identificar a cônica como uma elipse de semieixos 10 e 
2. Veja a Figura 31.1. 


Figura 31.1: A elipse G +57 =1. 


. Como a equação x? — 16y? + 8x + 128y — 256 não contém termos 


em xy, ela já se encontra diagonalizada, restando apenas comple- 
tar os quadrados em x e y: 


(xº + 8x) — 16(y? — 8y) — 256 
e EO [(y—-4)2-16]) =256 
(x+4)2— 16(y— 4)” = 16 
(x+4)2 (9-4)? 1 
6. ] 


Efetuando a translação 


( x =x+4 
y1 =y-—4, 


a equação que representa a cônica se transforma, no sistema de 
coordenadas x1y1, em 


Podemos identificar a hipérbole na Figura 31.2. 


Xv 
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>X 


2 


. ; end Foro 
Figura 31.2: A hipérbole | — + —=1. 


c. Reescrevendo a cônica 4x? — 20xy +25y? — 15x — 6y = 0 na forma 
matricial, obtemos 
vAv+Bv=o, 


+= (5 Je, 
) 


onde 


B=(-15 —6). 


A matriz À é a matriz simétrica que representa a forma quadrática 
g(x, y) = 4xº — 20xy + 25y? com respeito à base canônica. Não é 
difícil ver — os cálculos ficam para você — que os autovalores da 
matriz A são AM =0 e Ap = 29,e os autovetores normalizados 
são 


sds 


autovetor associado ao autovalor Ay = 0, e 
2/29 ) 


| ( Rr ) autovetor associado ao autovalor À = 29. 


Como (uy, u>! forma uma base ortonormal de R?, então 


Pot vis (SIND 2/VD 
=lm vol= (55/55 5/1455.) 
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é uma matriz ortogonal que diagonaliza a matriz A, e a matriz 
diagonal correspondente será 


00 
belo e 


com D = P'AP. Como det(P) = 1, a matriz ortogonal P repre- 
senta uma rotação em R?. 


Ee) 


e substituindo v = Pv, em v/Av, obtemos 


Considerando 


vVAv =— (Pvi)'A (Pv1) 
= vi (P'AP)vi 
=v'Dv,; onde P'AP =D 


“(3 2)(2) 


=00y%. 


A forma linear se transforma em 


Bv =B(Pvi) 
= BPvi 
= 5/29 2/00 Na 
=(5 =9 (5,05 545 )( on) 
=(-3VD o) ( 


= —3/29x1 ê 


X1 
)1 


Substituindo 
vVAv=29y e By=-3V29x 
em v'Av+By = 0, obtemos 
29y; —-329x, — 0. 
ou, ainda, 
NB, 


XM = 5 1 
1 341 
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onde identificamos facilmente a equação de uma parábola. Veja a 
Figura 31.3. 


XV 


Figura 31.3: A parábola x, = vB y2. 


5. a. Reescrevendo a equação 2xy +2xz+2yz— 6x — 6y — 42 = —9 na 
forma matricial, temos 


vAv+By=-9, 


onde 
x 011 
v=|y |, 4=|[1 01 eB=(-6 —-6 —4). 
Z 110 


A matriz A já foi diagonalizada no Exercício 3.b. Encontramos: 
1/46 
u=| —2/ V6 | autovetor associado ao autovalor À, = —1; 


1/6 
1/2 


= 0 autovetor associado ao autovalor À = —1: 


=1/v2 
1/v3 


us= | 1/y3 | autovetor associado ao autovalor À, = 2. 


1/3 
Como (u;, uz, uz) forma uma base ortonormal de Rê, temos que 
L/V6 1NZ 1/V3 


P=[u; Us us) = -2/V'6 0 1/3 
1//6 —1/V2 1/43 
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é uma matriz ortogonal que diagonaliza a matriz A, e a matriz 
diagonal correspondente será 


Vale também que D = P'AP. 


Observe que det(P) = 1, logo P representa uma rotação em Rº. 


Considerando 
X X1 
v=| y eps feios 
Z 2 


e substituindo v = Pv, em v/Av, obtemos 
vVAv = (Pvi) A (Pvi) 
= vi(P'AP)vi 


=viDvi : onde P'AP=D 


—1 00 X1 
= (x )1 Z1) O =" 0 y1 
0) 02 Z1 


= +24. 


Agora, substituindo v = Pv, em Bv, obtemos 
Bv = B(Pvi) 
= BPvi 


1//6 1/2 1/V3 x1 
—=(-6 -6-4) | —-2//6 0  1/V3 yi 
1/06 —1/V2 1/V3 z 


2 2 16 
="=X|——>)y1— —&=4. 
JC vo Na 
Portanto, substituindo 


2 à 
vAv=--y+2Z eBv=>x—>y—-—&a 
ERES O VD caia 
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na equação v'Avy + Bv = —9, obtemos 


== nda A EM 


2 2 16 
RT (a 


Completando os quadrados nas variáveis xy, Y1 € z1, obtemos a 
quádrica 


Leg) (aeb esfaedo) 


Agora, aplicando a translação 
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Ron 
y2 yr 
2 =2z7+ 


obtemos 
—5 +20 =1, 


que representa um hiperbolóide de duas folhas. 


. Reescrevendo a equação 2xy— 62x -+10/2y-+z—31 — 0 na forma 
matricial, temos 


vAv+Byv=31, 
onde 
x 010 
v=|[» |,4=[1 00 |[eB=(-6/2 10/21). 
Z 0 0 0 


Deixamos para você, novamente, o exercício de calcular os auto- 
valores e os autovetores correspondentes da matriz A. Obtemos: 


0 
u =| 0 autovetor associado ao autovalor Ay = O; 
1 
1/v2 
u=— 1 v2 autovetor associado ao autovalor À, = 1; 
0 
=] /N2 
us =— 1/ v2 autovetor associado ao autovalor À = —1. 
0 
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Como fuy, us, us) forma uma base ortonormal de Rº, temos que 


0 1/2 —1/v2 
P=[u; us, us] = 0 1/v2 1/v2 
1 0 0 


é uma matriz ortogonal que diagonaliza a matriz 4, e a matriz 
diagonal correspondente será 


00 
D=| 01 06 |; 
070 = 


onde D = PAP. 


Observe que det(P) = 1, logo P representa uma rotação em Rº, 
a saber, uma rotação de 7 K 4 radianos em torno do eixo-z. Con- 
siderando 


X X1 
V— y e vVi=— 1 ; 
Vá Z 


e substituindo v = Pv, em v/Av, obtemos 
vVAv = (Pvi)'A (Pv1) 
= vi (P'AP)vi 
— viDvi : pois PAP=D 


0 0 0 X1 
=(my72)/01 0 y1 

00 —1 Z 
=y—4- 


Agora, substituindo v = Pv, em Bv, obtemos 


Bv =B(Pvi) 
=— BPvi 


0 1/42 —1/42 x1 
=(-6V2 10V2 1) 0 1/42 1/v2 ] [ ] 
10 0 


z 
x1 

=(1 4 16) »1 ] 
za 


=x1+4y1/ + 1621. 


Portanto, substituindo 
VA,v =y-% e Bv=x1 +4y/+ 1671 
na equação v'Av + Bv = —9, obtemos 


vi—zi+x +4y + 167 =31. 


Completando os quadrados nas variáveis y, e z1, obtemos 
(vi +41) — (27 — 162) +m1 =31 


(1+22-4-(27—-8) +64 +27 =31 


e, por fim, 
x+29=-(y +22 +(27—8)2. 


Agora, aplicando a translação 


x =x1429 
»=)y1+2 
2=27—8, 


obtemos 
Lo o DR sd 
2=-»+a; 


que representa um parabolóide hiperbólico. 


.- Queremos montar a matriz A = (a;;), onde a;; = F(u;, u;), 
u =(1,0)eu=(1, 1). Temos: 


d11 =F(u,, u,) a AQE 0), (1, 0)) = 2: 
1 
da = F(u,, u,) =P 
( 


da = F(u,, U2) sl 


Logo, 


b. Queremos montar a matriz B = (b;;), onde b;; = F(vi, Vj), 
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vi=(2,1)evi=(1,-1). Temos: 
by =F(vi, vi) Elo To; (2, 1)) =3: 
bp =F(vi, v2) =F((2, D, (1,—1)) =9, 
ba = F(vo, vi) =F((1,-1), (2; 1)) = 0; 


bo2 =F(vo, v2) =F((1,-1), (1,—1)) = 6. 


Logo, 
3 9 
s-(5 E 


Expressando v| e v> em função de u; e uz (os detalhes ficam para 


2 


você), obtemos: 
vi =lu+Hl-u; 
6) =2-u, +(—1)-w, 


onde 
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Aula 3 2 


UM CAso PRÁTICO 


Ao final desta aula, você deverá ser capaz de: 


compreender os conceitos de autovalor e autovetor; 
reconhecer um escalar como autovalor de uma matriz; 
reconhecer um vetor como autovetor de uma matriz. 
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Um MODELO DE CRESCIMENTO 
POPULACIONAL 


Nesta última aula, vamos ilustrar como a teoria de autovalores e au- 
tovetores de matrizes com coeficientes reais pode ser usada para analisar 
um modelo de crescimento populacional. 


Iniciaremos nossa discussão com a apresentação de um modelo sim- 
ples de crescimento populacional. Para isso, vamos supor que certas 
espécies têm uma taxa de crescimento constante. Isso significa que a 
população cresce a percentuais iguais em intervalos de tempos iguais. 


Vamos considerar uma espécie em que cada indivíduo de uma gera- 
ção produz r novos descendentes e, logo em seguida, morre. Assim, se 
pn denota o número de indivíduos da população da n-ésima geração, 
supondo que as gerações se sucedem a intervalos de tempos iguais, 
temos que 


Pn = "Pn-1. 
Por exemplo, se r = 2, temos: po é a população inicial da espécie; 


pi =2po; 


De modo geral, temos p, = 2"po. E para r arbitrário, temos 
Pn = r"po. Esse modelo pode ser usado, por exemplo, para descrever 
a população de certa bactéria, na qual, a cada período de tempo, cada 
bactéria se divide em duas outras. Para esse modelo, a população cresce 
para o infinito se r > 1, decresce para zero se O < r < 1 e permanece 
constante se r = 1. 


Como você pode notar, esse modelo populacional é muito simples. 
Por exemplo, para a maioria das espécies o número de descendentes 
depende da idade dos pais. No caso da espécie humana, uma mulher 
com 50 anos de idade tem mais dificuldade de ter filhos que uma de 20 
anos. Estudaremos um modelo que leva em consideração esse tipo de 
complexidade. 


Vamos considerar uma certa espécie de pássaros em que o número 
de machos é igual ao número de fêmeas. Assim, basta controlar o 
número de fêmeas. Vamos supor, ainda, que o período de reprodução é 
de um ano e que, após o nascimento de uma nova fêmea, ela só poderá 


se reproduzir após um ano de vida. Antes de um ano ela será considera- 
da uma fêmea jovem e após um ano será considerada uma fêmea adulta. 
Podemos, então, denotar por: 


Pjn à população de fêmeas jovens após n anos (n períodos de repro- 
dução); 


Pa,n à população de fêmeas adultas após n anos. 


Vamos também assumir que, a cada ano, uma fração a de fêmeas 
jovens sobrevive e se torna fêmeas adultas, que cada fêmea adulta pro- 
duz k novas fêmeas jovens e que uma fração B de fêmeas adultas sobre- 
vive. 


A suposição de taxa de sobrevivência constante significa que a so- 
brevivência dos adultos independe da sua idade, o que nem sempre se 
aplica. 


Com as suposições anteriores, podemos relacionar a população de 
fêmeas jovens e adultas da seguinte forma: 


Pjn=KkpPjn-l 
Pan= GDjnA -BDan=1 


o que nos dá um sistema linear de ordem 2. Em notação matricial, 
podemos reescrevê-lo como 


P,=AP, 
onde 
n=(Go saco 6) 
Pan B 
Observe que 
P =APs; 


p=APR=A(AB)=A42Py; 
PB=AP=A(A?P) = AP); 
P=AP;=A(AÍPy) = A?Py, 


e, assim, de um modo geral, 


onde 
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é a matriz que representa a população inicial de fêmeas (jovens e adul- 
tas). 


Exemplo 32.1. 


Vamos considerar o modelo descrito anteriormente durante um pe- 
ríodo de 20 anos com matriz A dada por 


o 2 
42 (os 05) 


Essa matriz informa que cada fêmea adulta gera k = 2 fêmeas jovens 
a cada ano e que as taxas de sobrevivência são a = 0,3 para fêmeas 
jovens e B = 0,5 para fêmeas adultas. Observe que a < B significa 
que as fêmeas jovens têm menos chances de sobreviver que as adultas. 
Vamos supor, inicialmente, que temos 10 fêmeas adultas e nenhuma 


jovem; portanto, 
0 
n=(8,) 


Assim, após um ano, temos 


Dr SE 20 
Bean ron) lo) 


Como p;,; = 20 e pai =5, a população total de fêmeas é de 25 
indivíduos após um ano e a razão entre fêmeas jovens e adultas é 


Pjl 20 


4. 
Pal 5 


Após o segundo ano, temos 


0 2 20 10 
B=-48=(05 05 )( 5 )-(a5): 


O valor de 8,5 para fêmeas adultas pode ser interpretado como um 
total de 8 indivíduos. No entanto, como p;j2 = 10 e paz=8,5,a 
população total de fêmeas é de 18 indivíduos após dois anos, e a razão 
entre fêmeas jovens e adultas é 


Procedendo dessa forma, obtemos a seguinte tabela de valores: 


Tabela 32.1 


Fêmeas adultas Ea de fêmeas | Pjn / Pan 
Pa, n n+ Pa, n 


E 


E JR 

L| 20 [os | 2 | 
2/1 10 [| 8 [| 18º | 1] 
31 17 [1 70 71 /%2M0/| 238/| 
Ep E TOS MS Ee! 
[2 | 20% | 


E 
propos fps paso 
E E O E 


Retornando ao modelo geral, suponhamos que a matriz A tenha dois 
autovalores reais distintos, A, e À>, com autovetores correspondentes 
V| € v2, respectivamente. Como v; e v2 são linearmente independentes, 
eles formam uma base de R? e, portanto, podemos escrever 


Po=ajvi+Havo, com aj,az ER. 


Como B, = A” Py, temos que 


= Aº(ayvi +aovo) ; 


e, portanto, 


P,= ay4"vi +azA"vo. 


Agora, como v| é autovetor associado ao autovalor À 1, temos 
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Av, =Apvi; 


e, de um modo geral, Aºv, = Aj'v1. Analogamente, A”v> = Aj'v>. Por- 
tanto, podemos reescrever a equação 


P,=m4"vy +aA4"vo> 
na forma 
P=aApvi ta vo. 
aa ad ao O kN 
O polinômio característico da matriz A = o: é dado por 


p(x) = det(xhb — A) 
=x —Bx-ka, 


cujas raízes são 


A=5 (8440). 


Comok>0,0<a<1le0O<B<1I, temos que B? +40k > 0, 
portanto, a matriz A de fato possui dois autovalores reais distintos, À e 
À», como supusemos inicialmente. Vemos também que 


n=5(B+yBº+4oh) >0 


m=5(8-yB+4oh) <0, 
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ainda, que |A1| > |A2]. Assim, neste caso, o vetor P, pode ser reescrito 


como 
Pi= M avi + (x) arm : 
M 


Mo 
A 


seja, 2 x () quando n é muito grande. Nesse caso, teremos 


Agora, já que 


n 
< 1, temos que (2) —> 0 quando n > +oo, ou 


Passava. 


Isso significa que, após um tempo grande, a população fica propor- 
cional a vj. 


Exemplo 32.2. 


Dando continuidade ao Exemplo 32.1, como 4 = ( a = > 


temos que o polinômio característico é 


p(x) = -0,5x—0,6. 


Assim, os autovalores são 


di => (0,5+265) a 1,06 


me : (0,5 — 2,65) = — (0,56. 


Efetuando contas rotineiras que você pode conferir, obtemos os res- 
pectivos autovetores: 


O qa 1 
"= 0,53 E 


Observe, do autovetor v;, que 


l 
—— 21 
0,53 ea 
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o que explica a razão p;n/Pa,n na quinta coluna da tabela do Exemplo 
32.1. 


No exemplo anterior, trabalhamos com precisão de duas casas deci- 
mais nas aproximações numéricas. E claro que obteremos informações 
mais precisas se usarmos um número maior de casas decimais. 


Devemos, também, esclarecer algumas limitações desse modelo. As 
taxas de nascimento e morte de uma população de pássaros variam de 
ano para ano e, em particular, dependem do clima da região. Em nossa 
discussão, assumimos um meio ambiente constante. 


Muitos ecologistas também têm observado que as taxas de nasci- 
mento e morte variam com o tamanho da população. Em particular, a 
população não pode crescer mais depois de atingir um certo tamanho 
limite, pois incorre no problema da falta de alimento. E, ainda, se a 
população crescesse indefinidamente a uma taxa constante, ela iria su- 
perpovoar qualquer ecossistema. 


Exercício 32.1. 


1. Usando o modelo populacional desenvolvido neste capítulo, de- 
termine o número de fêmeas jovens e adultas após períodos de 1, 
2, 5, 10, 19 e 20 anos. Em cada caso, calcule também a razão 
Pjn/ Pan Considere 


m=( ni k=3, u=0,4e B=0,6. 


£3 Esperamos que você tenha apreciado os conhecimentos matemá- 
ticos desenvolvidos neste curso. Eles são, realmente, de ampla 
aplicação prática. Na medida em que você desenvolver outras 
ferramentas matemáticas, você verá esses conceitos ressurgindo 
em muitos contextos diferentes. No mais, nós, autores, desejamos 
a você toda a sorte e sucesso na sua caminhada pelo maravilhoso 
mundo da Matemática. 


TOTNAÇIN [E vanv 
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SOLUÇÕES DE EXERCÍCIOS SELECIONADOS 


AULA 19 


001 
ds no 1 ). 
100 


2. autovalor À, = 1 com multiplicidade 2: autovetores 
u, = (1//2,0,1/v2) eu, = (0,1,0); 
autovalor À) = —1 com multiplicidade 1: autovetor 


us = (1//2,0,-1/v2). 


AULA 20 


1. Matriz da projeção ortogonal com respeito à base canônica: 


a=( 1/4 3/4 E 


v3/4 3/4 
A diagonalização da matriz A é dada por 
A=PDP' = 
( 1/2 E or RA Ita Ee 
v3/2 1/2 0 0 —V3/2 1/2 
AULA 21 


Exercício 21.1 
0 0 
0 0 
0. 


1/05 0» 0 
x no 0 10 ] 
Los oro 


3. É dada pelo produto de matrizes 


da logs alo So di) 
A O q 
do do + E AUo a E FE 
E 01.0 5 ss l= 
o 000 La 2 

Va 46 vo vê Vê 
E CR E 
6 6 3 
1 5 1 

— 6 6 3 
“1 Lo 1 
3 3. 3 
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1. Como 4º = 4, temos 


(A?) = (AA) — AIM! Es (AM? =A2, 
garantindo que A? é uma matriz simétrica. 


2. Sejam P matriz ortogonal (P”! = P”) e D matriz diagonal tais que 
A = PDP". Então 


A? = AA = (PDP” (PDP') = PD(P'P)DP' = PDIDP' = PD?P, 


mostrando que A? também é diagonalizável por matriz ortogonal. 


3. Como A é uma matriz simétrica, temos, pelo Teorema 22.3, que 
A é diagonalizável por matriz ortogonal. Os autovalores de A são: 


A =3 com multiplicidade algébrica 2; 


Ap = —1 com multiplicidade algébrica 2. 


Uma base ortonormal para o autoespaço E (3) é dada por: 


u,=(1//2,1/2,0,0): 
= (0, 0, 1//2,-1/V2), 
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enquanto uma base para o autoespaço E(—1) é dada por: 


us = (1/42,0,0); 
w = (0,0,1/v2,1//2). 


Assim, as matrizes 


1/2. 0 1//2. 0 

E ou 1/2 0 —1/2 0 
0 1/v2 0 1/2 
O 1/2 0 1/42 


[a a a O) 
[a o SO a) 
| 
pa 


satisfazem 4 — PDP". 


AULA 23 


Exercício 23.1 


sto AGVDO MED No ao ça OO 
[; 1 P=-( a ha) o=(0 :) 
1/V3 —1/N2 —1/V6 
ne (hã 1/v2 ans) 
1/3 0 2/6 
= O 
»-[ 010 
001 
iva 0 ID O 
po. 1/02 0 1/02 0 | 
PEstEs O 1/2 o 1/2 |? 
O 1/2 0 142 
0000 
Our 0:00 
dad Ro 
0004 
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2. Observe que À = 5 é um autovalor de A, mas v = (—1,1,0) não é 
um autovetor correspondente ao autovalor À = 5. Temos: 


1/43 —1//2 —1//6 
P=[|1//3 1/v2 -1//6 |: D= 
1/3 0 2/v'6 


Sou 
SONO 
ND OO 


AULA ES MÓDULO 2 


AULA 24 


1. A matriz que representa o operador T com respeito à base canônica 


É 
Ze de d 
E DO 
| = 2 


Como A é uma matriz simétrica, segue que o operador T é au- 
toadjunto. 


2. A base pode ser 8 = (u,,u>,us|, dada por 


um =(—1/3,1/v3, 1/43); wo = (1/V2,1/V2,0) e 
U3 — (—1//6,1//6,-2/46). 


AULA 25 


1. Paratodou, v, weR”ecaeR, 
F(utaw,v) =(utaw)Av 
= (u +aw)Av 
=uAv+a(wAv) 
=F(u,v)taF(w,yv). 


Assim, F é linear na primeira variável. De forma análoga, mostra- 
se que F também é linear na segunda variável. 


L0/0 2 
ade es | O TO bB=|[ 1 
001 2 
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Exercício 27.1 


2 


: G % 
1. A hipérbole de equação Do qe l. 


AULA 28 


1. 2) =x5 —y5: parabolóide hiperbólico. 
2. x5 +y5 — 225 = — 19; hiperbolóide de duas folhas. 


2 2 2 
3. Z+2+52= 1: elipsóide. 


AULA 29 
Exercício 29.1 
LL Mh=2+ v=(=1+i,1) 
rh) =2-i) v= (—1 = 1) 
2. O polinômio característico é p(x) =x? —-2ax+ a? + b?, cujas 
raízes são À = a+bi e À = a — bi, com autovetores associa- 
dosvi = (1,-i)ev> = (1, i), respectivamente. 


3. Basta observar que, se 4 é matriz real, então seu polinômio ca- 
racterístico p(x) tem coeficientes reais. Logo, se À é uma raiz 
complexa de p(x), então À também é raiz de p(x). 
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1. Os autovalores são À = 1,44 e À = — 0,836, com autovalores 
correspondentes 
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Valores: 


Tabela 32.2 


Fêmeas jovens | Fêmeas adultas e de fêmeas | Pjn / Pan 
Pan n+ ER] 


a EE SS 
[EE DE CER 
Ea dE SE SE 


re ps fps pao 


o [600 
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